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INTRODUCTION 

The  experience  of  the  Analysis  Group  engaged  in  theoretical  si:j)ersonic  aero¬ 
dynamic  analysis  in  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins  Uni¬ 
versity  has  shown  that  the  computations  of  pressure  distributions  on  airfoils  of 
infinite  span  in  a  two-dimensional  flow  and  of  the  effects  of  shock  fronts  can  lx; 
greatly  speeded  up  by  the  use  of  appropriate  tables.  Brief  tables  for  the  calcula¬ 
tion  of  pressure  distributions  have  been  given  by  Lighthill1  but  so  far  as  was  known 
to  the  Analysis  Group,  when  its  work  was  started,  tables  to  facilitate  the  computa¬ 
tion  of  oblique  shock  front  effects  did  not  exist.  The  Analysis  Group  has  calculated 
extensive  tables  for  facilitating  oblique  shock  front  computations. 

At  the  same  time  the  need  for  a  concise  and  accurate  discussion  of  the  principles 
of  supersonic  aerodynamics  became  apparent.  With  the  aid  of  its  consultant.  Dr. 
F.  D.  Mumaghan,  Professor  of  Mathematics  in  The  Johns  Hopkins  University, 
the  Analysis  Group  prepared  such  a  discussion  for  two-dimensional  supersonic 
flows,  giving  careful  attention  to  the  precise  definition  of  terms.  An  important 
feature  of  this  study  is  the  demonstration  that  an  arbitrary  curve  intersecting  once 
and  only  once  each  of  the  straight-line  patching  or  characteristic  curves  in  an 
expansive  supersonic  flow  can  be  treated  as  if  it  were  a  characteristic  curve.  This 
fact  leads  at  once  to  the  formulas  for  calculating  pressure  distributions  over  airfoils 
immersed  in  two-dimensional  supersonic  flows,  and  hence  to  lift,  drag,  and  moment 
computations.  This  explicit  demonstration  is  not  given  in  any  of  the  literature 
that  so  far  has  come  to  the  attention  of  the  Analysis  Group.  Also  it  is  explicitly 
demonstrated  that  the  Ackeret  linear  theory  and  the  Busemann  second  order 
theory  are  approximations  to  the  rhore  general  theory. 

The  present  Bumblebee  report  includes  both  of  the  studies  just  described — the 
computation  tables  and  the  discussion  of  the  principles  of  two-dimensional  super¬ 
sonic  flow.  Thus  it  is  essentially  a  manual,  for  the  two-dimensional  case,  on  the 
theory  and  practice  of  computation  in  supersonic,  aerodynamic  flow,  as  far  as  the 
present  state  of  the  subject  is  concerned. 

1  F.  AI.  Panel  No.  G54,  Ac.  2)20;  A.  It.  C.  No.  7384. 
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PART  I.  GENERAL  THEORY 


1.  The  theory  of  'patching  curves  (characteris¬ 
tics)  in  two-dimensional  supersonic  flow. 

We  shall  confine  our  discussion  io  the  case  of  steady 
two-dimensional  flow  of  a  perfect  fluid.  The  basic 
equations  of  motion  are,  in  the  al>sencc  of  mass 
forces. 


1  3ji  _  du  1  Op  dv 

P  Ox  <it  ’  p  3v  (it 

We  |)ostulutc  a  relation 

!»  =  i»(p,  v) 

where  ij  is  the  entropy  and  assume  that,  rj  docs  not 
vary  with  x  and  y;  then 

3|>  3p  3 p  "  Op  3p  3p  3 p  *  Op 

Ox  Op  Ox  Os  ’  Oy  Op  dv  3v 

where 


so  that  e  is  the  “local  sound-velocity”.  Since  the 
flow  is  steady 


du 

dt 


du  du 

u  —  +  V  —  , 
Ox  av 


dv  3v  3v 
—  =  u  —  +  v  — 
dt  Ox  Oy 


and  so  the  equations  of  motion  may  lie  written  in 
the  form 


c-  Op  3u  3u 

- =  u - (■  v —  , 

P  Ox  Ox  dv 


c1  3p  3v  3v 

—  —  =  u  —  +  v  —  . 
p  dv  Ox  Oy 


The  equation  of  continuity  (for  steady  flow)  is 

3  .  ,  3  . 

—  (pu)  +  —  (Pv)  =  0. 

Ox  Oy 


Upon  multiplying  this  by  —  -  and  substituting  the 

P 

values  jus'  given  for  —  —  —  ,  —  —  —  we  obtain 

P  ox  p  dy 

(On  3u\  /  Ov  0v\  /  3u  0v\ 

"  s + v  r, ) +  X  s + v  »y ) -  c\s + *■) - 0 

i.e. 


(1) 


We  now  consider  an  arbitrarily  given  curve  C: 
x  =  x(a);  v  =  v  (a) 

(a  l)eing  any  convenient  parameter  or  independent 
variable)  and  suppose  that  the  values  of  u,  v  and  of 

the  vorticity  ^  are  prcscrilied  along  this 

curve.  In  other  words,  u,  v,  and  ~  arc  f?vcn 

functions  of  a;  does  it,  then,  follow  that  the  values 

of  the  four  space  derivatives  — ,  —  ,  v- ,  —  of  the 
*  .  dx  dy  dx  dy 

two  velocity  components  (u,  v)  arc  unambiguously 

determined  along  G?  It  not  we  can  regard  t  he  curve 

C  as  a  -patching  curve  along  which  two  different  flows 

are  patched  together,  there  l>cing  no  discontinuity 

in  the  velocity  components  (u,  v)  nor  in  the  vorticity 

t-  —  —  across  the  patching  curve.  For  this  reason 
Ox  dy 

any  such  patching  curve  C  must  !>e  clearly  dis¬ 
tinguished  from  a  shock-curve  across  which  the  nor¬ 
mal  component  of  velocity  is  discontinuous. 

Since  u  and  v  are  given  functions  of  a  along  C  we 
have 


3u  dx  3u  dy  du 

Ox  (hr  3y  da  da 

3v  dx  Ov  dy  dv 

Ox  da  Oy  da  da 


These  two  equations  together  with  the  equation 


Ov  0u 
Ox  Oy 


u(a), 


where  w(a)  is  a  given  function  of  a,  and  the  equation 
(1)  constitute  a  system  of  four  linear  equations  for 
the  four  space  derivatives  of  the  velocity  components 
(u,  v).  These  equations  determine,  without  am¬ 
biguity,  the  four  space  derivatives  of  u  and  v  unless 
the  determinant  of  the  coefficients  happens  to  be 
zero.  Thus  in  order  that  C  may  be  a  curve  which 
patches  together  two  flows  for  which,  at  the  various 

points  of  C,  at  least  one  of  the  four  derivatives  — , 

ox 


t-  ,  ~ ,  —  has  different  values  for  the  two  flows  (the 
o y  dx  dv 


values  of  u,  v  and 

dx 


—  being,  however,  the  same 

dy 


1 


i 
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for  both  Hous  along  C'«,  it  is  necessary  that 

This  quadratic  equation  for  the  ratio  4^  *  ^  has 

da  da 

real  roots  if,  and  only  if, 


+  v2  >  c2 

(the  roots  lwing  equal  when  the  equality  holds). 
Thus  a  point  of  the  flow  cannot  be  a  point  of  a  patch¬ 
ing  curve  unless  the  flow  is  supersonic  (or,  at  least, 
sonic)  at  that  point;  in  other  words  a  necessary 
condition  for  a  point  of  the  flow  to  be  a  point  of  a 
patching  curve  is 

q  >  c 

(where  q  =  (u2  +  v2)}  is  the  sj>ecd  of  the  flow). 
The  essential  difference  between  -supersonic  and 
subsonic  flow  is: 

There  are  no  patching  curves  in  subsonic  Jloic.~ 

It  remains  now  to  .see  what  further  condition,  if 

anv,  must  be  imposed  on  u,  v  and  —  —  —  (which 

dx  ch¬ 
are  given  functions  of  a  along  C)  in  order  that  the 

,  ..  ..  ..  du  du  dv  dv 

four  linear  equations  connecting  — .  —  mav 

dx  dy  dx"  dy 

have  a  solution  when  the  determinant  of  their  co¬ 
efficients  is  zero.  This  condition  is  obtained  by 
equating  to  zero  the  determinant  obtained  by 
replacing  any  one  of  the  four  columns  of  the  determi¬ 
nant  of  the  coefficients  of  the  four  equations  by  the 
column  of  numbers  on  the  right-hand  side  of  the  four 
equations.  Replacing  the  second  column,  for  ex¬ 
ample,  we  obtain 


dx  du 
da  da 

dv  dx 

0  r  ~ 

UU  UCK 


0 


w(a)  1 


0 

dy 

da 

0 


-  0 


c2  —  u2  0  —  uv  c2  --  v!l 

or,  equivalently. 

,  .  dx  uv  ,  ,  dy  du 

(c-  -  v2)  —  -  f  (c2  -  u2)  f  - 
da  da  da  da 

dx 


=  <  (c2  -  v2)  -f  uv 
da 


dy^  dx 
daj  da 


wfa). 


2  What  arc  here  termed  patching  curves  arc  usually  called 
characteristics;  \vc  prefer  the  name  patching  curve  since 
there  is  danger  of  confusing  characteristics  with  shock 
curves. 


We  shall  coniine  ourselves  to  the  case  where  «(a)  = 
0  so  that  the  vortieity  is  taken  to  be  zero  along  the 
patching  curve;  in  particular,  this  will  certainly  be 
the  case  when  the  flow  is  irrolalional,  i.e.  when  the 
vortieity  is  zero  throughout  the  entire  flow.  Then 
the  sufficient  condition  for  the  existence  of  solutions 
of  the  four  linear  equations  for  the  four  space  deriva¬ 
tives  of  u  and  v  along  C  (when  the  determinant  of 
these  four  equations  is  zero)  is 


dx  dv  dy  du 

(c2  -  v2)  —  —  -r  (c2  -  u2)  —  —  =  0. 
(let  (la  (la  da 


(3) 


From  (his  point  on  we  shall  suppose  that  the  flow 

is  irrolalional,  i.e.,  that  —  =  —  and  that  the  rela¬ 
te  dy 

tion  which  expresses  p  as  a  function  of  the  density 
p  and  the  entropy  v  is  of  the  form 

P  =  kpT 


where  k  is  a  function  of  the  entropy  tj  alone  (so  that 
k  does  not  vary  with  x  and  v).  Since,  when  the 
flow  is  irrotational, 


du 

dt 


3u  3v  3 

u - v  —  =  — 

Ox  Ox  Ox 


Uq;) 


dv 

dt 


3v 

u  —  +  v 
3x 


3v 

3y 


3u 

u  —  +v 
3v 


3v 

0y 


d_ 

dy 


(Rtf 


the  two  equations  of  motion  are  equivalent  to  (i.e. 
imply  and  are  implied  by)  the  single  relation 

r  c2  . . 

Jq2+  /  —  dp  =  const.  (4) 

J  P 


Using  the  relation  p  =  kpT  we  have  c2  =  ■—  =  ykp 


7-1 


=  —  SO  thll 


dp 


.  fc  .1  Tk  ,-i  T  p  C 
ip  7-1  7  —  1  P  7—1 

We  may,  then,  write  the  single  relation  which  is 
equivalent  to  the  two  equations  of  motion  in  the 
form 


C!  =  1  (q*j  -  qJ)  (5) 

where  qM  is  a  constant  ^  q  (it  is  the  value  of  q  at 
those  points,  if  any,  of  the  (low  where  c  =  0). 

The  differential  equation  which  determines  the 
patching  curves  may  be  written  in  the  foim 

(c2  —  v2)  (dx)2  +  2  uv  dxdy  +  (e2  —  u*)(dy)*  =  0 


0 
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and,  on  dividing  through  by  ds2  =  (dx)2  -j-  (dv)3, 
\vc  obtain 


On  denoting  by  8  the  inclination  of  the  flow  we  have 
u  =  q  cos  9,  v  =  q  sin  9 

so  that 

c1  /  .  dx  dv\: 

-  =  I  Sill  9- - cos 9~  I  -  Sill"  ft 

<1*  \  ds  ds  J 

where  y  is  the  angle  between  the  direction  of  flow 
and  the  direction  of  the  patching  curve;  y  is  the 
Mach  angle  and 

sin  it  —  c/q. 

The  number  M  =  q/c  =  cosec  y  is  known  as  the 
Mach  number  of  the  flow  at  the  point  in  question. 

Through  any  point,  then,  where  the  flow  is  sujrer- 
sonic  there  pass  two  possible  patching  curves  and 
the  direction  of  flow  bisects  the  angle  between  these 
two  curves  (the  angle  l>ctween  the  direction  of  flow 
and  the  patching  curves  1  icing  ±  y.  it  is  convenient 
to  direct  the  patching  curves  so  that  the  direction 
of  flow  bisects  internally  the  angle  between  the  two 
curves).  At  a  point  where  the  flow  is  sonic,  y  = 
<r/2  and  the  two  possible  patching  curves  through 
this  point  touch  at  the  point  (the  direction  of  flow 
being  along  their  common  normal  at  the  point). 

The  two  slopes  of  the  possible  patching  curves 
through  a  point  where  the  flow  is  supersonic  arc 

—  uv  ±  c\/q-  —  c- 
c'-  —  ii- 


Let  us  consider  the  curve  for  which 
dy  _  —  (uv  +  c\/ q-  —  c2) 


dx 


u- 


dq 


and  calculate  the  rate,  at  which  the  speed  q 


varies  along  this  patching  curve.  Since  q*  =  u'  +  v* 
and  since  (c"  —  v2)  dx  dv  +  (c2  —  u2)dy  du  =  0 
we  have 


q  dq  =  u  du  +  v  dv 


(  e2  —  u2dy 

=  <  u  —  v - - 


du 


=  <  u  +  v 


c-  —  v2dx 
uv  +  c  \/q"  -  c! 


du 


— ~  — -  (uc  +  V  v q*  -  c2)  du 


so  that 


(uc  -r  v  Vq1  —  c1) 


du 

da 


Since  u  =  q  cos  6,  v  =  q  sin  6,  we  have 


u  dv  —  v  du  =  ufsia  d  dq  -j-  q  cos  fi  dfi) 

—  v(cos  fi  dq  —  q  sin  fi  dfi)  =  q*dfi 

and  since 

.  ,  f  c1  —  u5dv  1  , 

udv  —  vdu  =  —  <  u -  —  —  v  r  du 

•  c~  —  v:dx  j 

=  { ~ — -  (uv  -f  c \/q*  —  e2)  -  vl 
\c-  -  V-  J 


dll 


v(q-  —  c2)  -f-  uc  \/ q-  — 


du 


c - v- 


we  obtain 


.dfi  Vq2  —  c*  - -  du 

q-—  =  - : - T  luc  +  v  Vq1  -  c1}  y  - 

da  c v*  da 

On  combining  this  with  the  formula  furnishing 

<l<t  . 
q  ^  we  find 


1  dq _ c 

q  dfi  V q2  —  C2 


tan  ft. 


(0) 


This  is  the  fundamental  relation  as  far  as  the  appli¬ 
cation  of  the  theory  of  patching  curves  (character¬ 
istics)  to  supersonic  calculations  is  concerned.  It 
is  easy  to  show  that  if  we  had  selected  the  other 
characteristic  direction  through  the  given  point  we 
should  have  obtained 


1  dq 
qdfi 


—tan  ft. 


Summarizing  these  results,  we  may  say: 

Along  any  patching  curve  (or  characteristic)  in  super¬ 
sonic  flow 


1  dq 
q  do 


=fc  tan  ft 


where  q  is  the  speed,  8  is  the  inclination  of  flow  and 
n  is  llic  Mach  angle. 

If  we  denote  by  t  the  component  of  velocity  tan¬ 
gent  to  the  patching  curve,  so  that  t  =  q  cos  y,  we 
have,  since  c  =  q  sin  y, 


Cl2  =  C2  +  i*. 

When  this  is  substituted  in  the  relation 


3 
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;  7  —  l  ,  ; 

c  =  2 —  Wa  ~  (D  we  gel 


T  +  1  .  7  —  1  ,  i 

__  ct  =  — — -  (q„  _  J-) 


so  that 


so  that 


c*  =  (fi“  ~  fI)- 

7  +  * 


Now  the  relation  t  =  q  cos  n  yields 

•It  =  cos  it  dq  —  q  sin  it  <!/• 

and  on  sulistituting  for  dq  its  value  it  q  tan  y  (16, 
we  obtain 


•It  =  q  sin  it  (±d B  —  «!>•)  =  c(±dP  —  dji). 


Hcnce 

±9 


-»+/7-«+ 4/^/7 iN 

j/5±Lm"  i 

V  7-i  <lu 


=  M  + 


+  const. 


Since 


sin-1  —  =  cot-* 

q* 


-  4/777- 

Thus  the  inclination  of  the  flow  at  any  point  of  a 
patching  curve  is  connected  with  the  Mach  angle 
by  the  formula 

±0  =  n  —  a  / tan-1  (  a  / 1  tan  #•  W  constant 

y T-i  \ty-i  / 

=  constant  —  fOi) 
where 

f to  =  j/: I±i  tan-'  ^  y^±I  tan  -  n .  (7) 

It  follows  from  this  relation  that  if  0  is  known  at 
every  point  of  a  patching  curve  and  if  y  is  known 
at  one  point  of  this  curve  then  y  may  be  determined 
at  every  point  of  the  curve.  All  that  is  needed  is 
a  table  of  values  of  the  function  f(jt).  For 

±(0  ~  6o)  =  }{n  o)  —  f(it) 


r(ft)  =  =b(0o  —  3)  +  f(jto). 


(S) 


This  relation  determines  f(y)  and,  hence  (by  means 
of  the  table),  y,  for  any  given  value  of  6. 

Table  5  furnishes  the  values  of  f(y),  calculated 
for  7  =  l.-J. 

7  -  1 


The  relation  c'  = 


(qji  —  q*)  implies  the 


relation  ({'  -  c5  =  7— ■  (<f  —  qf„)  where  q„,  = 


’7  -  1 


T  7  +  1 

is  supersonic, 


qM .  Hence,  at  any  point  where  the  flow 


q» 


>q  >  <i»  =  <!u  • 


In  a  velocity  diagram  (where  the  coordinates  of  any 
|x>int  arc  the  velocity  components  (u,  v))  a  point 
wlierc  the  flow  is  supersonic  lies  in  the  circular  ring 

q«  >  q  >  q« 

(the  polar  coordinates  of  any  point  in  the  velocity 
diagram  Ixiing  q  Z0).  The  relation 

1  dq 

-  —  =  db  tail  /» 
qd  3 

shows  that  the  graph,  in  the  velocity  diagram,  of  the 
points  on  a  patching  curve  lie  on  an  epicycloid  traced 
by  a  point  on  the  circumference  of  a  circle  of  diam¬ 
eter  qM  —  qm  which  rolls  externally  on  the  circle  of 
radius  q,„  whose  center  is  at  the  origin.  In  fact 
if  a  circle  of  radius  b  rolls  externally  on  a  circle  of 
radius  a  the  coordinates  of  a  point  on  its  circumfer¬ 
ence  are  given  by 


X  =  (a  +  1>)  cos-  t  —  l)  cosf  1  + 
»  \ 


y  =  (a  +-  b)  sin  -  r  —  b  sin {  1  H —  It 
a  \  a/ 

where  r  is  the  angle  through  which  the  circle  of 
radius  bjias  rolled.  The  polar  coordinates  r  Z  0  of 
this  point  are  given  by 

r2  =  x2  +  y2  =  (a1  +  2nb  +  2b2) 

—  2b(a  +  b)  cos  r,  tun  0  =  y/x. 

From  the  latter  relation  we  deduce  sec‘0  d0  = 

y  dx 


x  dy  -  y  dx  .  .  x  dy 

-  so  that  d0  = 


x- 


xs  +  i’¬ 


ll  dice 
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i4 1\9  =  x  ilv  —  v  dx 


(:i  -f  b)b  fr,  h  /  b\  *1 

!=  - - -  x  (:i  +  b)  COS-  -  —  h  cos^l  -r  -  Jr J 

r  h  (  b\  *1 

X  I  COS  -  r  —  cost  I  +  -  It  I 
+  £(:»  +  b)  sin  -  r  —  b  sin  ^  1  -j — ^ rj 
X  £sin  “  T  —  sin  ^  I  -f  ”^Tjl  *^T 


(a  +  b)b(:i  +  2I») 


(1  —  cos  t)  dr. 


In  the  case  of  the  velocity  diagram,  q  takes  the  place 
of  r  and  for  points  on  any  member  of  one  of  the  two 
families  of  patching  curves  we  have 


1  dq 

±-—  =  lan  it 
qdff 


_  - kfiH  Ji iLziil 

Vq1  -  c5  y  T  +  1  y  q*  -  qs.' 


Then  qM ,  the  radius  of  the  outer  circle,  takes  the 
place  of  a  +  2b;  and  q,„ ,  the  radius  of  the  inner 
circle,  takes  the  place  of  a.  We  calculate,  accord¬ 
ingly,  (a  -f  2b)2  —  r2  and  r2  —  a2: 

(a  +  2b)2  —  r2  =  2b(a  +  b)(l  +  cos  r) 
r2  —  a2  =  2b(a  +  b)(l  —  cos  r) 


so  that 


Since 


(a  +  2I>) 2  -  r2  _ 


r-  —  :i= 


T 

=  cot-. 
2 


r  dr  =  b(a  +  b)  sin  r  dr 


we  find,  on  using  the  value  of  r2dfl  calculated  above, 
that 


1  dr  _  a 
r  d0  a  +  2b 


cotI=  ;i .  /(a  +  2b)=  —  r= 

2  a  +  2b  y  r:  _  a* 

.  &  Qm  /y  —  1  .  r  ..  , 

Since  — r-rr  =  —  =  A/  f  t  it  follows  that  the 
a  +  2b  qM  y  7  +  1 

plot.,  in  the  velocity  diagram,  of  any  member  of  one 

of  the  two  families  of  patching  curves  (namely,  the 

family  for  which  ^  ^  =  +  tan  y)  is  an  epicycloid 

traced  by  a  point,  on  the  circumference  of  a  circle 
which  rolls  between  two  circles,  of  radii  q,„  and  qM  , 
both  circles  having  a  common  center  (the  origin). 


The  members  of  the  second  family  of  characteristics 
also  plot  as  epicycloids  traced  by  a  point  on  the  cir¬ 
cumference  of  the  same  rolling  circle:  all  we  have  to 
do  is  to  replace  r  by  —  r  (or,  equivalently,  by 
2x  —  -).  This  merely  amounts  to  the  geometric  fact 
that,  through  any  j>oint  in  a  circular  ring  there  can 
1)C  drawn  two  circles  each  of  which  touches  the 
circles  which  determine  the  ring  (the  common 
ratlins  of  the  constructed  circles  being  the  difference 
of  the  radii  of  these  two  circles). 

When  a  patching  curve  serves  to  patch  a  uniform 
flow  to  a  second  flow  (uniform  or  non-uniform)  the 
plot,  in  the  velocity  diagram,  of  the  patching  curve 
is  a  point  and  the  patching  curve  itself  is  a  straight 
line.  In  fact  u  and  v  are  constant  along  the  patch¬ 
ing  curve  so  that  the  plot,  in  the  velocity  diagram, 
of  the  patching  curve  is  a  ]>oint.  Since  q  is  constant 
along  the  patching  curve  so  also  is  c,  in  view  of  the 


relation  c2 


(q2  —  q2),  and  it  follows  from 


(2)  that  is  constant  along  the  patching  curve: 

in  other  words  the  patching  curve  is  a  straight  line. 
This  result  is  of  fundamental  importance: 

A  curve  which  yalchcs.  ill  mijwrsonic  flow,  a  uniform 
flow  onto  a  second  flow  {uniform  or  non-uniform)  is  a 
slraiijhl  line  whose  plot,  in  the  velocity  diagram,  is  a 
point. 

What,  can  we  say  about  the  various  memljcrs  of 
the  other  family  of  patching  curves  which  intersect, 
such  a  straight  line  patching  curve  (whose  plot  in 
the  velocity  diagram  is  the  point  P0  =  (u0,  v0))? 
The  plot  of  any  such  patching  curve  (of  the  second 
family)  in  the  velocity  diagram  must,  pass  through 
Po  and  hence  must,  be  the  epicycloid  (of  the  second 
type)  which  passes  through  P0 .  In  other  words: 

All  patching  curves  of  the  second  family  which  inter¬ 
sect  a  curve  which  patches  a  uniform  flow  to  a  second 
flow  {uniform  or  non-uniform)  have  the  same  plot  in 
the  velocity  diagram. 

Let  us  now  examine  the  various  members  of  the 
first,  family  of  patching  curves  which  intersect  these 
members  of  the  second  family.  The  various  points 
where  a  fixed  member  of  the  first  family  intersects 
different  members  of  the  second  family  all  have  the 
same  plot  in  the  velocity  diagram  (for  the  various 
members  of  the  second  family  plot  into  a  single  curve 
in  the  velocity  diagram).  lienee  the  velocity  vector 
is  constant  along  the  fixed  member  of  the  first  family 
of  patching  curves  and  this  implies  that  this  member 
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of  the  first  family  of  patching  curves  is  a  straight 
line.  In  other  words: : 

The  patching  curve  which  patches  a  uniform  flow 
to  a  second  flow  ( uniform  or  non-uniform )  is  not.  only 
a  straight  line  but  it  is  one  of  a  family  of  straighLlincs 
each  of  which  is  a  patching  curve;  furthermore  the 
velocity  vector  is  constant  along  each  member  of  this 
family. of  straight  lines. 

Along  each  member  of  the  family  of  straight  line 
patching '  curves,  both  q  and  6  sire  constant  so  t  hat 

the  relation  -  =  ±  tan  p  is  meaningless  (the  left- 

q  dO 

hand  side  of  this  relation  lx:ing  indeterminate).  If 

the  relation  -  ^  =  +  tan  u  is  valid  along  the  mcm- 
q  d0 

bers  of  the  second  family,  q  decreases  with  6  and 

hence  c,  p  and  n  increase  as  9  decreases.  If  -  -jr  = 
’ r  q  d0 

—  tan  p,  q  increases  as  6  decreases,  c,  p  and  p  de¬ 
crease  as  6  decreases,  and  the  flow  is  defined  as  an 
expansive  flow.  It  will  lx:  assumed  from  now  on 
in  this  paper  that  the  flow  is  expansive. 

Since  the  velocity  vector  is  constant  along  each 
member  of  the  family  of  straight-line  patching  curves 
so  also  is  c,  by  virtue  of  the  relation 

c  =  —TV—  (flM  -  q)- 


What  is  the  practical  significance  of  these 
facts?  If  we  draw  an  arbitrary  curve  (not 
merely  a  member  of  the  second  family  of  patching 
curves)  which  intersects  once  and  only  once  each 
member  of  the  family  of  straight-line  patching 
curves,  the  values  of  q,  0  and  p  at  any  point  of  this 
curve  are  the  same  as  those  at  the  corresponding 
point  of  any  member  of  the  second  family  of  patch¬ 
ing  curves  (two  points  being  said  to  correspond  when 
they  lie  on  the  same  straight-line  patching  curve). 
In  other  words  the  plot  in  the  velocity  diagram  of 
the  arbitrary  curve  is  the  same  as  the  plot  of  any 
one  of  the  members  of  the  second  family  of  patching 
curves.  This  implies  that  the  relation 
1  dq 


q  d(? 


=  —Ian  p, 


which  is  valid  for  any  member  of  the  second  family 
of  patching  curves,  is  also  valid  for  any  curve  which 
intersects  once  and  only  once  each  member  of  the 
family  of  straight-line  patching  curves.  In  par¬ 
ticular,  it  is  valid,  in  two  dimensional  flow,  for  the 
contour  of  an  obstacle  immersed  in  the  flow  if  each 


component  curve  of  the  contour  intersects  once  and 
only  once  each  member  of  the  family  of  straight-line 
patching  curves.  It  is  understood  that  8  does  not 
increase  as  we  move  from  the  nose  towards  the  trail¬ 
ing  edge  of  the  obstacle. 

To  avoid  two  sets  of  formulas,  one  set  for  the 
upix'r  part  of  the  contour  and  another  for  the  lower 
part  of  the  contour,  we  must  keep  in  mind  that,  as 
usual,  6  is  measured  positively  in  the  counter¬ 
clockwise  direction.  We  then  treat  the  lower  part 
of  the  contour  as  if  it  were  the  upper  part  of  a  con¬ 
tour  obtained  by  reflecting  the  given  contour  in  a 
line  through  the  nose  and  parallel  to  the  direction 
of  flow.  This  amounts  to  taking  9  positive  in  the 
clockwise  direction  when  we  are  working  on  the 
lower  contour.  We  now  have  a  basis  for  the  follow¬ 
ing  theorem  on  which  rest,  the  lift,  drag  and  moment 
calculations  in  siqxirsonic  flow: 

If  Pj  and  P  are  two  points  which  arc  both  on  the 
upper  part,  or  both  on  (he  lower  pari,  of  the  contour 
of  the  obstacle  then 

0,  —  0  =  fU)  —  Up) 


where  6,  p  are  the  inclination  and  the  Mach  angle, 


respectively,  at  P,  and  0i  ,  p\  are  the  inclination  and 
Mach  angle,  respectively,  at  Pi .  Since  0  and  0\  are 
given  by  the  geometry  of  the  obstacle  (they  are  the 
inclinations  of  the  tangents  at  P  and  Pi ,  respec¬ 
tively)  this  relation  furnishes  p  when  p\  is  given. 

Once  p  has  been  determined  it  is  easy  to  calculate 
the  ratio  p/pi  of  the  pressures  at  the  points  P  and 

Pi .  In  fact  the  relation  p  =  k py  yields  c"  =  -r-  — 

9p 

ykp7-1  so  that  -n  —  ( -\  and  —  =  ( -'j  = 

ci  W  Pi  W 

y 

i-' .  The  relations 

<y  —  |  ^ 

c  =  <i  sin  p;  e2  =  — —  (<1m  -  <15) 


yield 


q2(?  -  cos  2p)  =  {y  -  l)qj, 
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so  that 


C-  = 


4.  - 


2  I’  y  —  COS  2/i 


y  -  1  \  - 

><IM 


=  Cr-l)*K.(  ,-Y 

\7  —  COS  2m  / 


On  denoting  by  g(/i)  the  function 
follows  that 


\T  -  cos  2m/ 


in  gG<i) 

Table  5  furnishes  values  of  g(*x)  (calculate*!  for 
y  =  1.4). 

The  determination  of  the  pressure  at  any  point 
of  an  obstacle  in  supersonic  flow  is,  then,  merely  a 
matter  of  reading  tables  once  mi  is  known.  It  is 
convenient  to  take  Pi  as  the  nose  of  the  obstacle; 
once  mi  is  known  at  the  nose  f (m)  is  calculated  from 

fC«)  =  f(M.)  -  (0.  -  0) 


and  then  u  is  determined.  —  is  then  read  from  a 

Pi 

table.  Since  the  uniform  flow  usually  suffers  a 
shock  at  the  nose  of  the  obstacle,  mi  is  not,  in  general, 
the  Mach  angle  mo  of  the  undisturbed  flow.  Rather 
Mi  must  be  determined  by  a  separate  shock  calcula¬ 
tion,  the  theory  of  which  is  explained  in  the  following 
section. 

2.  The  theory  of  shock  surfaces  in  supersonic 
flow. 

The  conditions  which  hold  at  a  shock  surface  may 
readily  be  found  by  considering  the  time  derivative 
of  a  volume  integral 


I  =  f  f(x,  y,  z,  t)dV 
JV 

The  surface  S  bounding  the  volume  V  of  integration 
may  vary  with  t  and 


dl 

dt 


[  dV  +  f  fwn  dS 


where  wn  is  the  normal  component  of  velocity  of  the 
surface'  (the  normal  being  drawn  away  from  the 
volume  of  integration).  This  result  is  valid  if  f 
possesses  continuous  derivatives  with  respect  to 
x,  y,  z  and  t  throughout  V.  If  f  does  not  satisfy 


these  conditions  the  theorem  is  not  necessarily  valid. 
However  if  V  is  divided  by  a  surface  2J  into  two 
regions  Vo  thd  Vi  in  each  of  which  the  conditions 
hold  we  may  apply  the  theorem  to  V0  and  Vt  sepa¬ 
rately  and  combine  the  results;  the  only  penally  is 


that  we  have  to  consider  the  surface  integrals  of  fwn 
over  the  two  sides  of  We  obtain,  then, 

-  =  [  -dV+  f  fwndS  +  [  (f0  -  f,)w„  dr 

(the  normal  on  2  being  drawn  from  V0  to  Vj).  Ap¬ 
plying  this  theorem  to  a  volume  where  S0  and  Si  are 


composed  of  particles  of  the  fluid  and  are  such  that 
Vo  and  Vt  are  infinitesimal  we  obtain,  up  to  infinitesi¬ 
mals, 


dl 

dt 


I  fv„  dS  +  f 


fv„  dS  + 


fi)\vn  dr 


where  v  denotes  fluid  velocity  (as  opposed  to  w 
which  denotes  velocity  of  the  surface  2);  as  So  and 
Si  — >  2  we  obtain 


J  fo(w  -  Vo)„  dr  -  j  f,(\v  -  v,)n  dr. 

If  we  know  that  (y;  =  0  (for  every  V  whose  boundary 
(it 

consists  of  given  fluid  particles),  it  follows  (in  view 
of  the  arbitrary  choice  of  the  part  of  £  over  which 
the  integration  is  carried  out)  that 

fo('V  -  Vo) n  =  fl(\V  -  V,)„. 
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If,  in  particular,  2  is  stationary  we  find 

=  (f|V,)„. 

In  other  words  this  means  that  the  product  of  f  by 
the  normal  component  vn  of  the  fluid  velocity  Ls  con¬ 
tinuous  across  2.  This  is  certainly  the  case  (in  view 
of  the  principle  of  conservation  of  mass)  when  f  is 
the  density  p.  This  furnishes  the  firs!  shock  condi¬ 
tion  (for  stationary  shock  fronts),  namely 

(p»v«)a  =  (p,Vi)»  ,  (10) 

and  it  merely  expresses  the  principle  of  conservation 
of  mass. 

We  obtain  the  second  and  third  shock  conditions 
from  the  mechanical  law  that  force  =  time  rate  of 
change  of  momentum.  If  l  is  the  x-direction  cosine 
of  the  unit  outward-drawn  normal  we  have 

//'dS*i/v'ttdV' 

Setting  f  =  pu  and  proceeding  as  before  we  obtain 
J  (po  -  PiKdS  =  J  (poU«(wn  -  vn)0 

-  Piu,(wn  -  va),l  dS 

(the  outward  drawn  normal  to  S0 ,  for  example,  ap¬ 
proaches,  as  S  — »  2,  the  normal  to  2  drawn  from 
Va  towards  Vn).  Hence  at  any  point  of  2 

(po  -  Pi)£  =  poUo(wn  —  v„)o  -  P|Ul(\V„  -  v„),. 

Taking  the  x-axis  tangential  to  2  we  have  l  =  0  so 
that 

poUo(\Vn  —  Vn)o  =  PlUi(\V„  —  Vn)l  . 

On  combining  this  with  the  first  shock  condition  we 
obtain 

(vt)o  “  (vt)i-  (11) 

This  is  the  second  shock  condition;  the  tangential 
component  of  fluid  velocity  is  continuous  across  the 
shock  (whether  or  not  the  shock  is  stationary). 
Note  that  this  is  a  result  of  the  combination  of  two 
physical  laws:  a)  conservation  of  mass,  and  b)  con¬ 
servation  of  momentum  under  zero  force. 

The  third  shock  condition  follows  on  taking  the  x- 
axis  normal  to  2  so  that  l  —  1.  For  a  stationary 
shock  we  obtain 

Po  +  Po(vn)o  =  Pi  +  Pl(vn)g  .  (12) 


The  fourth  and  final  shock-condition  follows  from 
the  energy  principle.  The  rate  at  which  work  is 

being  done  cn  the  fluid  is  —  /  pv„dS  and  we  equate 

J* 

this  to 

ijbiw+unv 

where  U  is  the  internal  energy. 

This  gives 

Po(v„)o  +  p.(Jq;  +  U.)(v«  -  w»)0 

=  Pi(v«)t  +  pj(iqi  +  Ui)  (v,  -  w„), . 

For  stationary  shocks  this  reduces,  in  view  of  the 
first  shock  condition,  to 

>e  +  i<lo  =  ii  +  Iq? 

where  i  =  U  +  -  is  the  heal  content,  or  enthalpy. 

If  the  flow  is  isentropic  in  each  region  di  =  dU  + 

pd  +  —  =  —  and  writing  p  =  kpT  we  have 
\P/  P  P 

2 

di  =  yk p1-s  dp  so  that  i  =  — -  =  -  ?  (the 

y -  1 p  7-1 

constant  of  integration  being  determined  by  the 
convention  that  i  vanishes  with  p).  Now  the  argu¬ 
ment  which  showed  that  in  irrotational  two-dimen¬ 
sional  flow'  the  two  equations  of  motion  aie  equiva¬ 
lent  to  the  single  relation 

c5  =  y—^  (fl«  -  <ls) 

may  be  repeated  without  change  for  three-dimen¬ 
sional  flow;  we  find  that,  when  the  flow  is  irrotational, 
the  three  equations  of  motion  are  equivalent  to  the 
single  relation 

=  ~  (q>!  q5)  • 

Hence  the  shock  condition, 

+  ho  -  +  IQ?  (13) 

7-1  7-1 

may  be  formulated  as  follows: 

In  irrotational  flow  the  value  of  the  constant  qM 
is  the  same  on  both  sides  of  a  shock  surface.  I.n  two- 
dimensional  irrotational  flow,  then,  the  circular  ring 
in  the  velocity  diagram  whose  bounding  circles  have 
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radii  q#  and  qm  respectively,  is  unaffected  by  the  pres¬ 
ence  of  shock  lines  (straight  or  curved). 

In  view  of  the  relations  c2  =  — * ,  (v,)0  =  (vt)i 

P 

we  may  write  the  fourth  shock  condition  in  the  form 

But  the  first  and  third  shock  conditions  yield 

Pi  —  po  =  1  ~ ^ (vj;  =  pi 


so  that 

Ml  - 


.  p,  —  p0/ PI  po\  (pi  —  Po)  (pi  +  fif) 
o  (vJ  i  =  I  1 

Pi  po  \P0  Pi/  popi 

--(i  +-)--( i  +eY 

Pi  \  Po /  Po\  Pi/ 


Hence 

Pi/ 7+1  _  Pi\  _  Po/ 7+1  _  A>\ 
pi  \7  1  PV  Po\7  —  1  P>/ 

or,  equivalently, 

Pi  _  (7  +  i)pi  —  (7  -  Dpq  =  (7  +  1)*  —  (7  —  l) 

Po  (7  +  1)po  —  (7  —  Dpi  (7  +  1)  —  (7  —  l)x 

where  x  =  —  is  the  density  ratio  across  the  shock 

Po 

surface.  This  is  the  Rankinc-Hugoniot  relation.  It 
may  be  written  in  the  equivalent  form 

pi  =  (7+  Dpi+  (7  -  Dpj  =  (7  +  Dy  +  (7  -  D  .  . 

po  (7  —  Dpi  +  (7  +  Dpo  (7  —  Dy  +  (7  +  D 

where  y  =  -1  is  the  pressure  ratio  across  the  shock. 
Po 

Taking  p0  as  given  (7*  0),  x  =  —  cannot  surpass 

Po 

the  ratio  =  G  (with  7  =  1.4)  nor  lie  less  than 
7  -  1 

7  —  1  _  1 

7+1  G 

3.  Floio  against  the  nose  of  a  wedge. 

We  denote  by  w  <  ~  the  angle  from  the  direction 

£t 

of  the  incident  flow  to  the  upper  surface  of  the  wedge 
and  assume  a  plane  shock  inclined  at  angle  a  as 


shown.  Then  (v,)0  =  q0  cos  a;  (v„)0  =  q0  sin  a; 
(v,)i  =  qt  cos  (a  -  w):  (vn),  =  q,  sin  («  -  w). 
The  first  and  second  shock  conditions  yield 

poq*  sin  a  =  pifp  sin(«  —  w)  (10) 

qo  cos  a  =  »ii  cos  (or  —  w)  (17) 


Eliminating  qj  and  q0  by  division  we  obtain 

x  tan  (or  —  w)  =  tan  a  (IS) 

where  x  =  -  is  the  condensation  (so  that  x  lies  in 
Po 

the  range  -  <  x  <  6). 

The  relation  (18)  expresses  either  of  the  two 
variables  x  and  a  as  a  function  of  the  other  (the  angle 
w  being  supposed  given).  In  order  to  determine 
either  of  the  variables  x  and  a  (and,  hence,  the  other) 
we  need  a  second  relation  between  them.  This  is 
forthcoming  when  the  Mach  number  Mo  of  the  in¬ 
cident  flow  is  given.  In  fact  since  (v„)0  =  q0  sin  a, 
(v„)i  =  qi  sin  (a  —  w)  we  obtain  from  (12)  and  (1G) 


Pi  -  Po  =  Poq0  Sin 

Hence 


;in*a^l  —  =  p;qj  sin:  (a  —  w) (x  —  1). 


—  =  1  +  ^  sin'af  1  —  =  1  +  7M*  sin*  a(  1  — 

po  Po  \  x/  \  x/ 

since  c2  =  =  7"  (on  the  assumption  that  p  = 

dp  p 

kpv).  From  the  Rankine-Hugoniot  relation  (14)  we 
derive 

p.  _ _ 2T(x  -  1) 

Po  (7  +  D  -  (7  -  Dx 


and  so 


MnsinJa  = 


(7  +  1)  -  (7  -  Dx ' 
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Similarly 

M>  SM,S  (“  -  W)  =  (7+l)Xl  (y  ^T)  *  (20) 

In  order  to  facilitate  the  determination  of  the  ratio 
p,/po  when  Mo  and  w  are  given.  Tables  1-4  have 
been  prepared.  Their  use  is  discussed  in  detail  in 
Sections  8  and  10. 

Jf.  Approximate  evaluation  of  the  pressure 
ratio  p/pi. 

As  we  have  seen  in  the  previous  sections,  the  de¬ 
termination  of  the  pressure  ratio  p/po  involves  two 
separate  calculations.  In  the  first,  place  the  shock 
tables  furnish  the  pressure  ratio. pi/po  and  the  Mach 
angle  mi  after  the  shock  when  w  and  M0  are  given. 
Then  Table  5  furnishes  p/pi  at  any  j)oint  of  the 
contour.  On  multiplying  p/pi  by  pi/po  we  obtain 
the  desired  pressure  ratio  p/po .  The  shock  calcula¬ 
tion  is  the  more  troublesome  of  the  two  and  it  may  be 
completely  avoided  (if  w  is  not  too  large)  by  means 
of  the  following  observation.  It  may  be  verified3 
by  Tables  2  and  3  that,  when  w  is  not  too  large,  the 
relation 

f(/»i)  —  f(no)  +  W 

holds  to  a  high  degree  of  approximation.  Since 
w  —  9  =  f  Oil)  —  f(>0 

it  follows  that,  to  a  high  degree  of  approximation, 
-  e  -  fOo)  -  f GO. 

This  means  that  we  can  ignore  completely  the  shock 
calculation  and  evaluate  the  pressure  by  means  of 
Table  5  starting  from  any  point  P0  in  the  uniform 
flow  (the  value  of  6  at  P0  being  zero).  The  following 
samples  show  how  closely  this  “short-cut”  calcula¬ 
tion  approximates  the  calculation  made  by  using 
the  shock  tables.  It  was  convenient  to  base  this 
comparison  on  Mj ,  the  Mach  number  behind  the 
shock-wave.  The  result  of  the  “short-cut”  method 
is  marked  by  an  asterisk. 

3  For  example,  it  is  found  from  Table  1  that  Mo  =  1.777 
for  a  -  39°  and  w  =  5°,  so  that  mo  =  34.25  and  fGio)  = 
110.40  .  The  relation  in  question  yields  f (mi)  =  115.40°, 
whence  m  =  355.49°.  But  from  Table  2  it  is  found  that  when 
a  =  39°  and  w  =  5°,  Mi  =  1.605,  whence  mi  =  38.54°, 


M. 

=  2.13 

w 

-  1* 

2* 

5* 

10’ 

15*  20* 

22* 

24* 

25.07* 

M, 

-  2.001 

2.055 

1.316 

1.759 

1.5(73  1.340 

1.233 

1.096 

.9262 

M? 

-  2.091 

2.055 

1.946 

1.771 

1 .001  1.432 

1.361 

1.2S9 

1.250 

ilevia- 

tion 

-  0.0% 

0.0% 

0.0% 

0.7% 

2.4%  6.9% 

10.4% 

17.6% 

35.0%; 

M.  «= 

1.S5 

w 

*»  i* 

2* 

5* 

10* 

15* 

20* 

20.20* 

Mi 

=  1.SI5 

1.7S0 

1.676 

1.49S 

1.29S 

-9S02 

.917S 

Mi* 

=  1.S15 

1.7S0 

1.G7S 

1.509 

1.336 

1.140 

1.131 

ilevia- 

tion 

=  0.0% 

0.0% 

0.1% 

0-7%, 

2.9% 

16.3% 

23.2% 

Ml  1.40 

w  =  1* 

2* 

5* 

S* 

Ml  a  1.365 

1.330 

1.215 

1.075 

Mi*  =  1.365 

1.330 

1.217 

l.OSl 

rievia- 

tion  =*  0.0% 

0.0% 

0-1% 

0.6% 

The  behavior  of  the  ratio  Mf/Mi  shows  that  there 
is  a  high  degree  of  continuity  in  the  physical  variables 
as  w  changes  from  positive  through  zero  to  negative, 
since  the  simplified  method  then  becomes  exact. 
(When  w  is  negative  there  is  no  shock  and  the  flow 
is  expansive  from  the  very  beginning.)  It  should 
also  be  noticed  that  M?/Mi  seems  to  lie  almost 
independent  of  M0  for  wedge  angles  not  too  large. 

Two  further  properties  of  this  solution  are  worthy 
of  mention: 

(1)  The  simplified  method  gives  zero  entropy  in¬ 
crease  on  passing  through  the  shock  wave. 

(2)  The  calculations  using  the  shock  wave  tables 
cannot  in  general  give  a  description  of  conditions 
in  the  wake,  since  there  are  two  independent  physi¬ 
cal  quantities  to  be  matched  on  the  two  sides  of  the 
wake,  and  only  one  adjustable  degree  of  freedom, 
namely  the  direction  of  the  final  flow.  Using  the 
simplified  method,  the  difficulty  is  of  a  different 
sort  since  now  an  infinite  number  of  solutions  exist. 
Since  p,  p,  M  depend  only  on  M0  and  on  the  angle 
between  the  incident  and  final  flow  directions,  the 
two  flows  derived  from  the  upper  and  lower  surfaces 
always  match,  no  matter  what  direction  is  assumed 
for  the  final  uniform  flow. 

5.  The  linear  theory. 

It  is  easy  to  derive,  from  the  rather  complicated 
formulas  of  the  preceding  sections,  simpler  formulas 
which  serve  as  approximations  which  arc  valid  for 
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sufficiently  small  values  of  w  and  of  0.  If  w  is 
so  small  that  its  square  is  negligible,  we  have 

Ism  («  —  w)  =  tstfi  «  —  w  secret 

so  that,  from  formula  (18) 


sec1  a 

x  =  1  +  w - -  1  +  2a-  cosec  2a. 

ism  a 

Hence  as  w  — » 0,  x  — >•  1,  and  formula  (19)  then  shows 

that  sin1  a  — »  rjj  =  sin2  mo  •  Thus  a  first  approxi- 
iVla 

mation  to  the  shock  angle  a  is  the  Mach  angle  mo 
of  the  undisturl)cd  uniform  flow.  Hence  a  second 
approximation  to  x  is 

x  =  1  +  2w  c-oscc  2>i,  .  (21) 

It  follows  from  formula  (14)  that  a  second  approxima¬ 
tion  to  the  pressure  ratio  pi/p«  across  the  shock  line  is 

Pi  _  1  +  (7  +  l)w  coscc  2mo 

pc  1  —  ('  —  l)w  cosec  2>io 

Since  we  are  neglecting  squares  of  w  this  may  l>c 
written  in  tl,e  equivalent  form 

—  =  I  +2y  \v  coscc  2/io .  (22) 

Po 

It  follows  from  formula  (20)  that  as  w  — *  0, 
Mi  — >  Mo  so  that  a  first  approximation  to  mi  is  mo . 
In  order  to  obtain  a  second  approximation  to  mi 
we  must  first  obtain  from  formula  (19)  a  second 
approximation  to  the  shock  angle  a.  To  do  this 
we  set  a  —  nn  +  kw  (where  k  is  a  multiplier  which 
we  wish  to  determine)  and  x  =  1  +  2w  coscc  2p0 
in  (19).  Since  Mo  =  cosec  m»  we  obtain 


Thus  a  second  approximation  to  mi  is 


Mi  =  mj-t- 


7  —  cos  2/i<> 

-  ,A- 

2  cos*  /io 


Now  it  follows  from  formula  (7)  that 


f  '00  = 


2  COS1  M 
y  —  COS  2/i 


(24) 


(25) 


so  that  to  a  second  approximation 

f(/o)  =  !(mo)  +  w.  (20 

This  is  the  basis  for  the  short-cut  calculation  of 
Section  4. 

It  follows  from  formula  (25)  ami  the  relation 
6  —  w  =  f(p)  —  f(/ii)  that 


—  w  =  (/i  —  It i) 


2  cos1  mi 
y  —  cos  2/ii 


Hence 


m  —  Ml  = 


7  —  COS  2/1, 
2  COS1  Ml 


(e  -  W) 


7  —  COS  2/ip 
2  ccs*  /id 


(0  -  w) 


(27) 


since  mo  is  a  first  approximation  to  m  ■  On  com¬ 
bining  (24)  and  (27)  we  obtain 


M 


MO 


+ 


7  —  COS  2/m  „ 

- - ff 

2  COS*  Mo 


(28) 


On  taking  the  logarithmic  derivative  of  the  func- 
/  sin2  m  \  J2. 

lion  gW  =  (7  _  tos  j if-'  "cob,ara 

-  -2-Ll »  -  -^r)  -  2r  •  W 

g(/l)  y  -  1  (  7  —  COS  2/1 J  7  —  COS  2/1 

It,  follows  from  formula  (9)  that, 


1  +  2k\v  cot.  no  = 


1  +  2iv  coscc  2/ij 
1  —  (7  —  l)w  coscc  2 Mo 


o  Cl)l  Mi _ 

7  —  cos  2 mi 


(m  -  Ml) 


so  that  k  = 


=  1  +  (7  +  1)\V  coscc  2/10 
(7  +  1)  COSCC  2pn  _  T_+-l 


Hence 


2  cor,  Mo  -I  cos'  Mo 

a  second  approximation  to  the  shock  angle  a  is 


,  y  +  l 

o  =  mo  +  - — ; — w. 
4  cos-  /io 


(23) 


On  setting  mi  =  Mo  +  k'w  in  formula  (20),  after  re¬ 
placing  Mi  by  37—  and  making  use  of  (23)  and  (21) 


sin  mi 


we  find  that  k' 


7  +  1 

2  COS"  Mo 


cos  2mo 


2  cos2  mo 


and  this  reduces,  by  means  of  formula  (27),  to 

—  =  1  +  2y  coscc  2m*  •  (0  —  w) 

I>i  (30) 

=  1+27  COSCC  2 Mo-  (0  —  w) 

since  we  arc  neglecting  w(0  —  w).  On  combining 
formulas  (22)  and  (30)  we  obtain 

—  =  1  -j-  2y0  COSCC  2m  1) .  (31) 

Po 

This  simple  formula  for  the  pressure  ratio  —  is 

Po 
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known  as  the  linear  theory.  Since  e2  =  — *  fonr.ula 

p 

(30)  may  lxj  rewritten  in  the  following  equivalent, 
form: 

PoV  _  Co  Pa 

Ap  =  p  -  p0  =  - - e  =  - - : - 9 

sin  Mi  cos  mo  sm  pa  cos  pj 

9  .  6 

=  pa  (lorn  : - -  poq0  — —  W-) 

Sin  Mo  COS  Mo  COt  Mo 


—  Po  q5 


VM|  -  1 


6.  The  quadratic  or  Busemann  theory. 

This  theory  is  obtained  from  the  method  of  this 
report  by  carrying  out  the  various  approximations 
as  far  as  ihe  w2  terms.  Since 

tan  (a  —  \v)  =  tan  a  —  see1  a-  \v  4-  sec5  a  tan  a-w5 

we  have 

tan  (a  —  w)  see1  a 

- =  I - w  -f  sec5  a-w5 

tan  a  tan  a 

SO  that 


so  that 

— —  =  . — - r-~  I  (7  +  I)  cjs  2 m  cose c5  p  —  4  cos1  m)  ■ 

g(p)  (7  -  cos  2m)- 

I  Icnee 

p  .  .  2y  cot  Ml  ,  x  ,  _  y_ 

-=14 - —  'P  -  Pi)  +  , -  r. 

Pi  7  —  CDS  2pi  (7  —  cos  2 Mil- 

X  1(7+  1)  cos  2pi  coscc2  pi  ~  4  cos5  pil 

X  (m  —  Pi)1  *  •  ■  (35) 

On  differentiating  (25)  we  obtain 


f"(p)  = 


-2(7  +  1)  sin  2p 


1  w  ~  1  o 

(7  -  cos  2m)- 

and  so  the  relation  0  —  w  =  f(p)  —  f(pi)  yields 

2  cos5  Pi  ,  t  .  (7  4-  1)  sill  2p,  ,  x. 

g  w  =  - —  (p  —  Pi)  +  - r~—  tp  -  Ml)-- 

7  —  cos  2pi  (7  —  cos  2pi)- 


Invcrting  this  relation  we  obtain 

7  cos  2pi  v  .  (7  +  l. 

p  -  mi  =  (f)  -  w)  4- - 


2  ccs5  pi 


(7  -f  1)  sin  2pi  (7  —  cos  2pi) 
S  cos*  Ml 

X  (0  -  w)s 


tan  a  sec5  a  ( see'  a  \ 

X'=  - -. - -  =  I  4 - W  +  I - see-  a  I  W- 

tan  («-w)  tan«  \tan5a  /  and  on  substituting  this  in  (35)  we  find 


sec5  a 

=  1  4 - w  4-  coscc5  a-w5. 

tan  a 

On  substituting  for  a,  in  the  coefficient  of  w,  the 
second  approximation 

-y  *{■  1 

a  =  Po  +  : - 2 —  w  (see  (23))  we  obtain,  since 

4  cos  po 

y  ~b'  f 

tan  a  =  tan  uo  +  : — ri —  w 
4  COS  Mo 

sec5  a  1  sec5po  ,  7  4-  1  ,,  ..  N 

- =  tan  a  4- ; - -  — - - 1-  - - - —  (1  -  cot-  po)w. 

tan  a  tan  a  tan  mo  4  cos’  mo 

Hence  the  third  approximation  to  x  is 
x  =  1  4-— — °w 

<7  (33) 

4-  {  coscc5  po  4-  - -  (2  -  coscc5  mo)  r  w* 

[  4  cos'  po  J 

On  substituting  this  value  of  x  in  (14)  we  obtain 
^  =  1  4-  27  cosec  2p0-\v 

P°  (  /  x  i  (34> 

,  f  .  ,  4  4-  (7  -  3)  ccscc5  p0\  „ 

4-74  cosec5  po  4- - : - : - r  w*. 

(  4  cos'  po  j 

On  differentiating  (29)  with  respect  to  m  we  obtain 

g*(p)  _  /gTp))* 
g(p)  l  g(p) ) 


(7  —  cos  2 p)5 


1  (7  --  cos  2 p)  cosec5  p  4-  2  col.  p  sin  2p) 


n  ,  .  y  sec*  pi 

-  =  1  4-  27  coscc  2pi  ( 0  —  w)  4 - ; — 

Pi  4 

X  [(7-  4-  1)  coscc5 pi  —  4  cos5  pil (0  —  \v)5 .  (36) 

We  replace  in  this  formula  mi  by  p-0  in  the  coefficient 

of  (0  -  w)2  and  bv  +  7  C^-—  w  (see  formula 

2  cos-  po 

(24))  in  the  coefficient  of  0  -  w.  In  this  way  .we 
obtain 

-  =  1  4-  27  coscc  2p0-  (0  -  w)  -  27  scc5p0  {coscc  2p0  cot.2p0  X 
Pi 

(  7-  cos2p0)  }w(fl  —  w)  4-  -  sec' po| (7  4-  l)cosec5po 

-  4  cos5pol(0  -  w)5 

=  1  4-  27  cosec  2po-  (,6  —  w)  4-  7  scc5po(7  coscc5 2p0 
4-  27  coscc  2pocot.  2po  —  eot52po)w5 
—  27  scc5po(7  coscc5  2p0  4-  7  coscc  2pocot  2po)wO 
4-  7sec5po(7  coscc2 2p0  r  cot52po)05- 
Upon  multiplying  this  by  (34)  we  obtain 

-  =  1  4-  2-,-  9  coscc  2p0  4-  7  see2  po  coscc2  2p0  X 

P«  (37) 

(7  4-  cos22po)05- 

This  approximation  to  the  pressure  ratio  is  known 
as  the  quadratic  or  Busemann  Theory. 


PART  II.  APPLICATION  OF  THE  GENERAL  THEORY 


7.  Limitations  on  method. 

This  part  of  the  report  describes  in  detail  the 
calculation  of  the  pressure  at  any  point  of  an  airfoil 
in  a  supersonic  stream.  The  lift,  drag  and  moment 
may  then  lie  obtained  by  numerical  integration, 
when  the  pressure  is  known  at  enough  points  on  the 
surface  of  the  airfoil.  In  making  the  application 
it  should  be  remembered  that  the  method  is  subject 
to  the  following  limitations: 

(J)  The  effect  of  viscosity  in  forming  a  boundary 
layer  is  neglected.  It.  is  assumed  in  effect  that  the 
flow  past  the  surface  is  not  retarded  by  frictional 
forces. 

(2)  The  calculations  are  strictly  two-dimensional. 
No  attempt  is  made  here  to  correct  for  the  effect 
of  finite  span. 

(3)  The  procedure  breaks  down  for  wedge  angles 
which  are  too  large.  For  a  given  incident  Mach 
number,  there  is  a  maximum  wedge  angle  through 
which  the  flow  may  be  turned  by  a  plane  shock  wave 
issuing  from  the  corner.  For  greater  wedge  angles, 
the  shock  is  presumably  curved  and  detached  from 
the  corner.  As  applied  to  airfoil  calculations  this 
limitation  means  not  only  that  the  leading  edge  of 
the  airfoil  must  be  fairly  acute  (which  is  likely  to  be 
true  of  any  airfoil  worthy  of  consideration),  but  also 
that  the  angle  of  attack  must  not  be  too  great.  No 
matter  how  thin  the  airfoil,  this  limits  the  permissi- 
ble  range  in  angle  of  attack.  This  limitation  is 
more  restrictive  for  lower  Mach  numbers;  in  fact, 
it  is  always  possible  for  the  Mach  number  to  be  so 
low  (yet  greater  than  unity)  that  the  calculations 
cannot  be  carried  out  for  any  angle  of  attack.  As  an 
example,  for  the  symmetrical  bi-con  vex  GU2  profile 
with  a  semi-angle  of  11.42°  at  the  leading  edge,  the 
computations  can  only  be  carried  out  for  angles  of 
attack  up  to  13.43°  for  M0  =  2.13°,  and  up  to  8.45° 
for  M9  =  1.85°;  for  M0  less  than  1.4S9,  the  calcula¬ 
tions  cannot  be  carried  out  for  any  angle  of  attack 
including  zero. 

(4)  This  method  is  not  exact-  even  for  the  ideal 
two-dimensional  non-viscous  case.  Unless  the  ob¬ 
stacle  is  a  wedge  formed  by  two  half  planes,  the 
shock  wave  is  necessarily  curved  and  the  motion 
has  vorticity.  However,  the  initial  slope  of  the 


shock  wave  is  given  by  the  present  method,  and  the 
curvature  of  the  shock  wave  is  small  if  the  curvature 
of  the  airfoil  is  small;  hence,  the  effect  of  vorticity 
will  ojxjrate  chiefly  near  the  trailing  edge  and  will 
probably  l)e  small.  It  seems  probable  that  this 
effect  is  small  in  comparison  with  the  effect  of  the 
boundary  layer  near  the  trailing  edge,  which  is 
likewise  ignored. 

In  all  of  these  tables  the  ratio  of  the  specific  heats 
(7)  is  taken  as  1.4.  The  formulas  on  which  the 
tables  are  based  are  restated  lxdow  for  this  value 
of  7. 

S.  Plane  shock  wave  solution  for  flow  in  a 
corner. 

The  following  formulas  descrilx:  two-dimensional, 
non-viscous  supersonic  flow  in  a  corner,  throughout 
the  range  of  Mach  numbers  and  corner  angles  for 
which  the  shock  wave  is  plane.  By  juxtaposition 
of  two  such  corner  flows  with,  in  general,  different 
wedge  angles  we  .have  the  case  of  a  uniform  stream 
impinging  on  an  infinite  wedge,  not  necessarily 
placed  symmetrically  with  respect  to  the  stream. 
On  either  side  of  the  shock  wave  the  flow  is  uniform. 
The  angle  w  through  which  the  flow  must  turn  in 
order  to  stay  in  the  corner  is  called  the  wedge  angle. 
The  Mach  number,  pressure  density  and  velocity 
are  denoted  by  M,  p,  p,  q,  with  subscript  0  or  1 
used  to  refer  respectively  to  the  state  before  or  after 
the  shock  wave.  The  angle  which  the  shock  wave 
makes  with  the  incident  flow  is  denoted  by  a.  The 
relations  between  these  quantities  are  as  follows: 


PI  tail  a 
po  tan  (a  —  w) 


13 


Two-Dimensional  Supersonic  Pressure  Calculations 


Tables  i-4  give  x,  M0 ,  Mi ,  and  —  respectively  as 

IV)  - 

functions  of  a  and  \v.  The  usual  problem  is  to  find 

Mt  and  —  ,  given  w  and  Mu .  Since  \v  and  Mo  are 
IV) 

known,  a  may  lx;  found  from  the  table  of  M0  as  a 

function  of  a  and  \v  (Table  2).  Then  Mi  and  — 

Po 

may  Ik;  read  from  Tables  3  and  4.  In  general  a 
two  way  interpolation  is  required  in  each  table.  A 
sjxxiific  numerical  illustration  will  be  given  later. 

Usually  the  x-table  (Table  1)  is  not  used.  It  was 
a  necessary  step  in  the  preparation  of  Tables  2-4, 
and  is  included  here  because -it  gives  a  quantity  of 
physical  interest,  namely  the  ratio  of  densities 
behind  and  in  front  of  the  shock  wave. 

Several  properties  of  this  solution  may  lx;  pointed 
out. 

The  density  is  greater,  the  Mach  number  smaller, 
the  pressure  greater,  the  velocity  smaller,  and  the 
velocity  of  sound  greater  behind  the  shock  wave 
(subscript  1)  than  in  front  of  the  shock  wave  (sub¬ 
script  0). 

The  shock  angle  «  is  approximately  equal  to  the 
Mach  angle  p  for  small  wedge  angles,  as  would  lx; 
expected  on  general  grounds.  (In  general  small 
wedge  angles  give  weak  shocks.) 

For  a  given  Mach  number,  there  is  a  maximum 
wedge  angle  w„, .  For  wedge  angles  greater  than 
wm  the  problem  has  no  solution  in  terms  of  a  plane 
shock  wave;  presumably  the  shock  moves  forward 
from  the  edge  and  is  curved.  There  is  another, 
slightly  smaller,  critical  wedge  angle  w,  for  which  the 
Mach  number  Ms  is  just  equal  to  unity;  for- wedge 
angles  between  ws  and  wm  the  flow  behind  the  shock 
wave  is  subsonic.  In  this  case  the  procedure  of  the 
next  two  sections  is  not  applicable  since  the  Mach 
angle  is  not  real;  furthermore,  it  seems  likely  that 
the  solution  is  unstable  even  for  the  case  of  an  infinite 
wedge,  if  the  flow  behind  the  shock  wave  is  subsonic. 
It  therefore  seems  likely  that  the  critical  angle4  for 
which  the  shock  breaks  away  from  the  nose  is 
nearer  to  ws  than  to  w)n  .  Practically,  the  distinc¬ 
tion  between  wm  and  ws  is  quite  unimportant.  These 
angles  are  tabulated  in  Table  0.  They  approach 
equality  for  both  small  and  large  Mach  numbers; 

1  It  is  not  certain  that  such  a  critical  angle  exists.  Since 
the  plane  shock  wave  solution  is  not  really  “exact”,  the 
shock  may  he  detached  slightly  and  have  some  curvature 
even  for  small  wedge  angles. 


(he  maximum  difference  w„.  —  ws  is  0.420°  (occur¬ 
ring  for  Mo  =  1.52).  As  Mc  increases,  w„,  and  w„ 
both  increase;  but  it  should  ne  noted  that  for  any 
Mach  number,  no  matter  how  large,  w„,  and  w,  are 
each  less  than  45.584°. 

Table  (5  was  constructed  from  the  following 
formulas  given  by  CroccoJ. 


With  cither  subscript, 

w  -  a  -  0, 

tan  0  2  /  sin1  jio  ^  7  —  1 

tan  a  7  +  1  \  sin*  a  2 


An  alternative  formula  may  be  derived  for  the 
determination  of  wm  ,  by  eliminating  a  between  the 


first  two  of  equations  (38).  Writing  2  = 


m  =  Mo2 , 


az*  +  bz  +  c  =  0 

where  a  =  4(6m  +  5)3(m  +  5) 

b  =  27 (m  4-  5)2  —  (m  —  l)’(6m  +  5)J  + 
I8(m  +  5)(m  —  l)(6m  5) 

c  —  — 4(m  —  l)3 

A  similar  formula  may  be  obtained  for  ws ,  also  of  the 
fourth  degree  in  tan  w, ;  but  it  cannot  lx;  solved 
explicitly. 

There  is  a  wedge  angle  w*,  analogous  to  ws  or  w,„  , 
which  limits  the  existence  of  the  “short  cut  method” 
described.  Since  f(p)  =  f(juo)  +  w,  and  since 


Mo 

1 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1,7 

1.8 

1.9 

2 

w* 

0° 

1.3' 

3.8* 

8.2' 

9.0* 

11.9' 

14.9” 

17.8* 

20.7’ 

23.6’ 

26.4' 

W* 

c* 

1.4' 

3.7' 

8.3' 

9.0' 

11.7' 

14.2’ 

10.6’ 

18.8' 

20.9' 

22.7” 

Wm 

o'" 

1.5” 

3.9' 

8.7' 

9.4* 

12.1' 

14.6' 

17.0' 

19.2' 

21.2° 

23.0” 

2.5 

3 

4 

5 

6 

10 

00 

39.1' 

49.8” 

65.  S' 

70.9' 

S5.0° 

102.3' 

130.5' 

29.7° 

34.0' 

38.7' 

41.1' 

42.4° 

44.4" 

45.6' 

29.8'’ 

34.1' 

38.8' 

41.1' 

42.4° 

44.4' 

45.6’ 

5  L.  Crocco,  Singolarita  della  corrcnte  gassosa  ipcracus- 
tica  nell’intorno  di  una  prora  a  diedro.  I/aerotecnica. 
Vol.  XVII,  fasc.  6,  1937. 
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Application  of  Theory 


f(*x)  ^  90°  (\/o  —  l)  =  130.-15°,  \v  must  Ik;  loss 
than  130.45°  —  f(pn)  if  m  is  to  Ik*  real.  Therefore 
w*  =  130.45°  —  f(p„).  The  preceding  table  shows 
that  w*  is  of  the  same  order  of  magnitude  as  ws  and 
wm  for  small  Mach  numtars,  but  diverges  radically 
at.  very  high  Mach  numlxrrs. 

Inspection  of  Table  2  shows  that.,  for  fixed  w  and 
varying  a,  M«  decreases  as  a  increases  up  to  a  point 
in  the  neighliorhood  of  05°,  after  which  it.  increases. 
Thus  for  each  wedge  angle,  there  is  a  critical  value 
of  Mo ,  below  which  the  plane  shock  wave  solution 
does  not  exist./’  In  Table  2,  the  smallest  value  of 
Mo  for  each  wedge  angle  is  enclosed  in  a  box.  Simi¬ 
larly  the  value  of  M0  for  which  the  flow  behind  the 
shock  becomes  sonic  is  underlined  and  in  italics. 
For  any  wedge  angle  and  for  any  Mo  above  the 
minimum,  two  values  of  a  may  lx;  read  from  Table 
2,  or  at  least  could  lx;  if  the  table  were  extended. 
The  smaller  of  these  two  values  of  a  should  always 
lx;  chosen.  The  solution  obtained  on  using  the 
other  value  involves  subsonic  flow  lx;hind  the  shock 
wave  and  is  probably  unstable.  In  any  case,  the 
method  of  this  report  is  applicable  only  to  super¬ 
sonic  flow. 

9.  Flow  around  a  corner  or  convex  curve  (“ex¬ 
pansive  flow'’). 

The  Mach  angle  of  the  initial  uniform  flow  is 

Pi  =  sin-1  — .  The  Mach  angle  p  =  sin-1  ~ 

at  any  point  where  the  flow  Las  turned  through  an 
angle  r\  is  determined  by 

f(pi)  -  f(p)  =  >1 

where  f(p)  =  Vo  tan-1  (V 6  tan  p)  —  p. 

Table  5  gives  the  numerical  value  of  f(p)  in  de¬ 
grees.  Since  f(p)  is  a  monotonically  increasing 
function  of  p,  it  is  seen  that  p  steadily  decreases, 
and  hence  M  steadily  increases  as  we  move  with 
the  flow.  Since  p  is  an  essentially  positive  quan¬ 
tity,  the  maximum  angle  through  which  the  flow  may 
turn  is  rjmnx  =  f(pi).  Practically,  however,  this  is 
no  restriction  since  f(pi)  is  greater  than  a  light, 
angle  for  Mach  numbers  less  than  2.5,  and  greater 
than  ()0°  for  Much  numbers  less  than  4.3. 

6  The  relation  between  this  critical  Mo  and  w  is  precisely 
the  same  as  the  relation  between  Mo  and  wm  given  by  Table 
6.  In  constructing  Table  6,  a  very  slightly  different,  ap¬ 
proach  was  made,  in  that  the  maximum  w  was  found  for 
fixed  Mo  (instead  of  the.  minimum  Mo  for  fixed  w). 


After  p  has  lxxm  determined  at  any  point,  the 
pressure  (p),  velocity  (q)  and  density  (p)  are  de¬ 
termined  by  the  following  formulae: 


p_ 

pi 


P  _  gOO  ,  .  _  /  sin 'y.  Y 

Pi  ~  g(«i)  K  M  \0.2  +  sin -nj 


0.2  4-  sii^Pi 
0  2  -f  sinsp 


hW 

fiM 


where 


h(p)  =  lg(p)lsfl 


/  sin5p  \s,: 
\0.2  +  siirp  / 


It  was  convenient  to  use  the  above  definition  of 
g(p),  which  differs  by  a  (non-essential)  constant 
from  the  function  g(p)  defined  in  Section  1.  It  is 

to  lx;  noted  that  p,  p,  and  ^/jp  (tl,c  local  velocity 

of  sound)  decrease,  and  q  increases,  in  flow  around 
a  corner  or  convex  surface.  Of  the  functions  f(p), 
g(p),  h(p)  tabulated  in  Table  5,  only  f  and  g  arc  used 
in  obtaining  the  pressure.  The  function  h(p)  was 
a  necessary  step  in  calculating  g(p).  and  is  given 
in  the  table  since  it  may  occasionally  lx;  of  interest 
to  know  the  density  changes  in  expansive  flow. 


10.  Specific  numerical  illustration. 

The  calculations  for  two-dimensional  airfoils  are 
based  on  the  two  problems  just  discussed.  The  first 
problem,  with  the  aid  of  Tables  1—1,  gives  the  pres¬ 
sure  and  Mach  number  on  either  surface  of  the  airfoil 
at  the  leading  edge,  just  behind  the  shock  wave  (the 
pressure  and  Mach  number  in  t  he  undisturbed  stream 
being  given).  The  Mach  number  thus  obtained 

is  converted  into  a  Mach  angle  (y.  =  sin-1  — ^ 

and  the  flow  around  the  curved  profile  is  followed  by 
means  of  the  solution  to  the  second,  problem.  In 
this  way  it  is  possible  to  calculate  a  value  for  the 
pressure  and  the  Mach  angle  at  any  point  of  the 
profile.  The  shock  waves  proceeding  from  the 
trailing  edge,  and  the  conditions  in  the  wake,  arc 
ignored,  on  the  principle  that  they  do  not  affect  the 
pressure  on  the  airfoil. 

As  a  numerical  example,  consider  a  C1U2  profile 
(formed  by  two  22.84°  circular  arcs),  and  specifically 
the  30%  chord  station  on  the  upper  surface,  with 
+5°  angle  of  attack  and  Mach  number  1.77.  The 
semi-angle  at  the  leading  edge  is  1 1 .42°,  and  this 
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Two-Dimensional  -Superson  ic  Pressure  Calculations 


would  be  the  wedge  angle  to  use  in  connection  with  a 
zero  angle  of  attack;  but  as  a  result  of  the  5°  angle 
of  attack,  the  wedge  angle  to  l>e  used  for  calculations 
on  the  upper  surface  is  w  =  11.42°  —  5°  =  0.42°. 
The  shock  angle  a  is  then  to  be  found  from  Table  1 
by  double  interpolation.  Inspection  of  Table  2 
shows  that  M0  =  1 .77  and  w  =  0.42°  determines  a 
to  bo  between  40°  and  41°.  For  w  =  6.42°,  interpo¬ 
lation  gives  Mo  =  .58  X  1.779  +  .42  X  L822  = 
1.797  for  a  =  40°  and  M„  =  .58  X  1.740  +  .42  X 
1.7S3  =  1.758  for  a  =  41°.  Another  linear  inter¬ 
polation  shows  that  for  M0  =  1.77,  a  =  40  +  %%  = 
40.09°.  Now  that  a  is  known,  Mi  and  pi/p«  are 
obtained  by  interpolating  in  Tables 3  and  4. 

Mi  =  .31  X  (-5S  X  1.573  +  .42  X  1.5S0) 

+  .09  X  (.5S  X  1.536  +  .42  X  1.541)  =  1.550 

Pi/Po  =  .31  X  (.5S  X  1.359  +  .42  X  1  -434) 

+  69  X  (.5S  X  1.354  +  .42  X  1.429)  =  1.387 

This  describes  the  state  just  behind  the  shock  wave. 
To  follow  the  subsequent  changes,  find  the  Mach 
angle  /xi  corresponding  to  the  Mach  number  Mi: 

sin  in  =  M,  =  .6452 


in  =  40.18°, 

and  determine  n  at  any  point  by 
t(Mi)  —  (GO  =  V 

where  ij  is  the  angle  between  the  tangent  at  the  point 
in  question  and  the  tangent  at  the  leading  edge. 
For  example,  the  geometry  of  the  wing  shows  that 
at  the  30%  chord  station  jj  =  0.88°,  so  that  at  this 
point  (using  Table  5) 


f (m)  =  f(40.1S°)  -  6.8S°  =  117.0S0  -  6.8S°  =  110.20° 
M  =  34  +  -,W  =  34.10° 

g 00  -17843 

g(Ml) 


Then  ~ 


=  .7045.  This  is  the 


Pi  gUnj  .25328 
ratio  of  the  pressure  at  any  point  to  the  pressure  just 
behind  the  shock  wave.  To  find  the  ratio  of  the 
pressure  to  the  pressure  in  the  undisturbed  stream, 
this  must  be  multiplied  by  pi/p0 .  In  the  above 
example,  the  ratio  of  the  pressure  at  the  30%  chord 
station  to  the  pressure  in  the  undisturbed  stream  is 


p/po  =  1-387  X  .7045  =  .9771 


The  calculations  for  the  lower  surface  proceed 
similarly,  but  with  a  wedge  angle  w  =  .11 .42°  -f-  5°  = 
16.42°. 

After  the  pressure  has  been  obtained  at  a  sufficient 
number  of  points  on  both  upper  and  lower  surface, 
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the  lift,  drag,  and  moment,  forces  are  obtained  by 
numerical  integration.  For  the  GU2  profile  it  is 
sufficient  to  calculate  the  pressure  at  10%  chord 
intervals,  including  both  leading  and  trailing  edges, 
and  integrate  by  Simpson’s  rule.  For  profiles  in  the 
shape  of  a  polygon,  the  integrations  may  be  effected 
without  approximation.  In  any  case,  the  integrals 
to  be  evaluated  are 


J  p  tan  0  dx 
J  (x  +  y  tan  0)  p  tlx. 

on  both  surfaces.  Here  x  is  the  distance  back  along 
the  chord  (measured  from  any  fixed  point  on  the 
chord  line;  the  moment  derived  is  then  understood 
to  l>e  the  moment  about  this  fixed  origin),  y  the  dis¬ 
tance  of. the  point  on  the  surface  above  the  chord 
line,  j8  the  angle  which  the  tangent  to  the  surface 
makes  with  the  chord  line.  All  of  these  quantities 
must  lie  taken  with  correct  sign:  y  is  positive  on  the 
upper  surface  and  negative  on  the  lower  surface;  f)  is 
positive  on  the  front  of  the  upjier  surface  and  the 
rear  of  the  lower  surface,  but  negative  on  the  rear 
of  the  upper  surface  and  the  front  of  the  lower  sur¬ 
face;  x  takes  on  negative  values  unless  the  origin  is 
taken  at  the  leading  edge.  The  moment  per  unit 
span  is 


/  (x  +  y  tan  0)p  dx  -  (x  +  y  tan  /3)p  dx 

^uppcr  ♦'lower 

which  is  considered  as  positive  if  it  tends  to  rotate 
the  wing  about  the  origin  in  the  sense  of  increasing 
angle  of  attack.  The  lift  and  drag,  measured 
perpendicular  and  parallel  to  the  chord  line,  are 


X  =  /  p  dx  —  /  pd> 

■Mower  ♦'upper 

x  -  /  p 
•'ll  pper 


-  p  tan  /3  dx 

♦'lower 


tan  /3  dx 


per  unit  span.  The  lift  and  drag,  measured  perpen¬ 
dicular  and  parallel  to  the  undisturbed  stream,  are 
per  unit  span 

L.—  X  cos  <j>  —  S  sin  <#> 


D  =  X  sin  <f>  —  5  cos  <f> 
where  <p  is  the  angle  of  attack. 


Application  of  Theory 


Although  this  procedure  is  occasion-ally  referred  to 
as  an  “exact”  method,  it  is  actually  an  approxima¬ 
tion,  although  presumably  a  very  good  one.  Two 
facts  must  lx;  neglected  in  the  development  of  this 
method: 

(1)  The  mathematical  solution  for  the  flow  around 
a  corner,  or  around  a  cylindrical  surface  which  is 
never  concave  outward,  is  based  on  the  assumption 
that  the  characteristic  lines  proceed  from  the  lxkly 
to  infinity  without  meeting  any  barrier  or  locus  of 
discontinuity.  In  the  wedge  problem,  all  of  the 
characteristics  on  lx)th  sides  of  the  shock  wave  inter¬ 
sect  the  shock  wave.  Hence  the  flow  around  a 
corner  (or  convex  surface)  cannot  lx:  treated  inde¬ 
pendently  of  a  preceding  flow  in  a  corner  involving 
shock.  Likewise  the  shape  of  the  shock  wave  is  not 
independent  of  the  flow  following  the  shock  wave. 
The  theory  of  the  plane  shock  wave  is  developed  on 
the  assumption  that  the  magnitude  and  direction  of 
flow  just  behind  the  shock  wave  is  everywhere  the 
same.  This  is  impossible  if  expansive  flow  follows. 
It  appears  certain  that  for  an  airfoil  such  as  the  GU2, 
a  general  theory  taking  account  of  this  interaction 
would  lead  to  a  curved"  shock  wave,  perhaps  slightly 

7  There  arc  physical  reasons  for  Iwlicving  that,  at  dis¬ 
tances  which  arc  large  compared  with  the  maximum  diam¬ 
eter  of  a  body  disturbing  a  uniform  flow,  the  shock  angle 
approaches  the  Mach  angle  of  the  undisturlicd  stream. 
This  leads  to  a  shock  which  continually  decreases  in  in¬ 
tensity  with  distance  from  the  disturbance,  as  compared 
with  the  present,  solution  which  gives  a  shock  wave  which 
has  the  same  “strength”  at  every  point. 


detached  from  the  leading  edge,  even  for  small  wedge 
angles;  in  this  ease  the  critical  wedge  angle  discussed 
earlier  would  have  no  exact  meaning. 

(2)  The  conditions  in  the  wake  have  been  ignored, 
presumably  on  the  theory  that  the  wake  is  lxmnded 
by  shock  waves  proceeding  from  the  trailing  edge, 
and  hence  can  have  no  effect,  on  the  pressure  at  any 
point  on  the  surface.  In  general,  however,  this 
method  docs  not  yield  a  solution  for  the  flow  in  the 
wake,  since  it  is  necessary  to  match  three  physical 
quantities  such  as  p,  p,  M,  of  which  two  arc  indc- 

» t"' 

pendent  (lx:eause  of:  the  relation  —  - ;  = 

2  y  —  1 

constant.)  and  there  is  only  one  degree  of  freedom  in 
applying  the  method  of  this  report,  namely  the 
direction  of  the  final  flow.  Therefore  the  assump¬ 
tion,  that  conditions  in  the  wake  do  not  affect  the 
flow  over  the  airfoil,  is  without  any  satisfactoiy 
foundation. 

Neglecting  these  two  effects  might  lxv  justified, 
partly  by  photographic  evidence  and  partly  by  the 
consideration  that  the  pressure  distribution  is 
affected  mainly  near  the  trailing  edge,  where  the 
pressure  is  small,  and  therefore  the  total  lift,  drag 
and  moment  arc  affected  only  slightly.  However, 
both  of  these  faults  in  this  method  should  lx:  recog¬ 
nized  as  mathematical  inconsistencies,  quite  distinct 
from  objections  that  may  be  raised  on  physical 
grounds  (e.g.,  the  neglect  of  viscosity). 
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TABLE  ! 
tan  a 


a 


w 

* 


ft 


X 


tan  (a-w) 


1° 

2° 

3° 

4* 

5° 

8° 

7’ 

8° 

9° 

10° 

11° 

12° 

13° 

14° 

16° 

16° 

17° 

18° 

19° 

20°  v. 

\m 

1° 

2.001 

1.501 

1.334 

1.251 

1.201 

1.168 

1.145 

1.127 

1.113 

1.102 

1.094 

1.C86 

1.080 

1.075 

1.070 

1.066 

1.063 

1.060 

1.05? 

1  J1.063 

2° 

3.002 

2.002 

1.669 

1.603 

1.403 

1.337 

1.290 

1.256 

1.227 

1.206 

l.iee 

1.173 

,  1-161 

1.150 

1.141 

1.133 

1.126 

1.12) 

1  1.133 

3* 

4.008 

2.5C6 

2.  CCS 

1.7S3 

1.606 

1.507 

1.436 

1.383 

1.342 

1.309 

1.263 

1.261 

1.242 

1.226 

1.213 

1.201 

1.190 

1  ii.213 

4° 

5.012 

3.010 

2.343 

2.010 

1.810 

1.678 

1.583 

1.S12 

1.458 

1.414 

1.378 

1.349 

1.324 

1.303 

1.285 

1.269 

lfl.303 

3° 

; 

6.021 

3.516 

2.682 

2.266 

2.015 

1.549 

1.731 

1.643 

1.574 

1.520 

1.475 

1.438 

1.4C7 

1.381 

1.358  , 

V'?. 

1  1.407 

6° 

7.034 

4. 025 

3.022 

2.522 

2.222 

2.022 

l.eeo 

1.774 

1.692 

1.626 

1.573 

1.529 

1.491 

1.480 
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1.525 

1.563 

1.584 

1.020 

:  .  1.383 

1.387 

1.391 

1.396 

1.403 

1.409 

1.417 

1.426 

1.436 

1.446 

1.459 

1.472 

1.467 

1.504 

1.523 

1.544 

1.568 

1.596 

1.628 

1.662 

1.702 

1.432 

1.436 

1.441' 

-  1.447 

1.453 

1.461 

1.469 

1.479 

1.490 

1.502 

1.515 

1.530 

1.547 

1.565 

1.586 

1.610 

1.630 

1.667 

1.701 

1.740 

1.786 

;  1.483 

1.487 

1.492 

l;499 

1.506 

1.514 

1.523 

1.533 

1.545 

1.558 

1.573 

1.589 

1.607 

1.628 

1.651 

1.677 

1.706 

1.739 

1.777 

1.820 

1.870 

1.535 

1.540 

1.545 

1.552 

1.560 

1.568 

1.578 

1.589 

1.602 

1.816 

1.632 

1.660 

1.609 

1.692 

1.717 

1.746 

1.777 

1.813 

1.864 

1.902 

1.966 

1.590 

1.596 

1.600 

1.607 

1.615 

1.624 

1.635 

i.647 

1.600 

1.675 

1.693 

1.712 

1.733 

1.757 

1.7S4 

1.815 

1.849 

1.868 

1.933 

,1.984 

2.044 

;  1.647 

1.661 

1.657 

1.664 

1.673 

1.682 

1.693 

1.706 

1.720 

1.737 

1.755 

1.775 

1.798 

1.824 

1.853 

1.886 

1.923 

1.965 

2.013 

2.069 

2.1® 

{  1.705 

1.730 

1.716 

1.723 

1.732 

1.742 

1,754 

1.767 

1.782 

1.799 

1.819 

1.841 

1.806 

1.893 

1.924 

1.969 

1.9® 

2.044 

2.096 

2.155 

2.223 

I  1.767 

1.771 

1.777 

1.766 

1.794 

1.804 

1.816 

1.830 

1.646 

1.064 

1.885 

1.908 

1.934 

1.963 

1.996 

2.034 

2.070 

2.124 

2.179 

2.242 

2.316 

1.831 

1.835 

1.841 

1.048 

1.867 

1.868 

1.881 

1.895 

V.  912 

1.931 

1.962 

1.977 

2.004 

2.035 

2.070 

2.110 

2.155 

2.206 

2.206 

2.332 

2.410 

"j  1.898 

1.902 

1.907 

1.915 

1.924 

1.936 

1.948 

1.963 

1.960 

2.000 

2.022 

2.048 

2.077 

2.110 

2.147 

2.188 

2.236 

2.2® 

2.352 

2.424 

2.506 

>  1.967 

1.971 

1.976 

1.983 

1.992 

2.004 

2.017 

2.032 

2.050 

2.071 

2.094 

2.121 

2.152 

2.186 

2.225 

2.269 

2.319 

2.376 

2.442 

2.517 

2.®  5 

|  2.041 

2.043 

2.048 

2*055 

2.064 

2.075 

2.089 

2.105 

2.123 

2. 145 

2.189 

2.197 

2.229 

2.285 

2.305 

2.362 

2.405 

2.46S 

2.534 

2.613 

2.706 

I  2. 117 

2.119 

2.124 

2.130 

2.139 

2.130 

2.164 

2.  ISO 

2.199 

2.221 

2.2(6 

2.275 

2.308 

2.346 

2.388 

2.437 

2.492 

2.556 

2.626 

2.711 

2.eoe 

l  2.198 

2.199 

2.203 

2.359 

2.217 

2.228 

2.211 

2.258 

2.277 

2.300 

2.326 

2.366 

2.390 

2.429 

2.474 

2.525 

2.682 

2.6(9 

2.725 

2.812 

2.914 

2.283 

2.286 

2.286 

2.291 

2.899 

2.309 

2.323 

2.339 

2.369 

2.382 

2.409 

2.440 

2.475 

2.516 

2.662 

2.615 

2.676 

2.745 

2.824 

2.916 

3.022 

|  2.373 

2.371 

2.373 

2.377 

2.384 

2.394 

2.407 

2.424 

2.443 

2.467 

2.494 

2.526 

2.563 

2.605 

2.eS3 

2.708 

2.771 

2.843 

2.926 

3.022 

3.134 

I  2.467 

2.464 

2.464 

2.467 

2.474 

2.483 

2.498 

2.512 

2.532 

2.556 

■^584 

2.616 

2.654 

2.698 

2.747 

2.805 

2.870 

2.946 

3.032 

3.132 

I  2.568 

2.563 

2.561 

2.663 

2.568 

2.576 

2.589 

2.054 

2.624 

2.648 

2.677 

2.710 

2.749 

2.794 

2.846 

2.904 

2.972 

3.060 

3. 140 

3.244 

3.366 

I*  2.675 

2.667 

2.663 

2.663 

2.667 

2.675 

2.686 

2.701 

2.721 

2.745 

2.774 

2.806 

2.847 

2.893 

2.946 

3.007 

3.078 

3.159 

3.252 

3.3® 

3.4® 

k  2.788 

2.778 

2.772 

2.770 

,2.772 

2.778 

2.788 

2.803 

2.822 

2.846 

2. STS 

2.909 

2.950 

2.997 

3.06J 

3.114 

3.187 

3.271 

3.368 

3.4® 

3.611 

f;  2.910 

2.898 

2.887 

2.883 

2.883 

2.887 

2.898 

2.910 

2.928 

2.962 

2.901 

3.015 

3.056 

3.105 

3.161 

3.225 

3.300 

3.367 

3.487 

3.604 

3.740 

£  3-M0 

3.022 

3.010 

3.002 

3.000 

3.002 

3.010 

3.CE2 

3.040 

3.083 

3.091 

3.126 

3.166 

3.217 

3.274 

3.341 

3.418 

3.608 

3.611 

3.732 

3.873 

TABLE  2 

u  (  5X  1 

M0  =  v  g  _  v  j  2  cosec  o  ] 


a°  !  22*  !  23°  24°  SB°  I  zt* 


»*  29°  I  an® 


ae®  > 


4.086  3.779 
iJOt  3.960 
4.407  4.146 
4.717  4.369 


6.032  4.33E 
6.414  4.946 
*•6*7  6.382 
6.631  6.823 
7.426  6.476 


6.828 

11.46! 

20.748 


3.366  3.182  3.023  2.692  2.764  2.637  2.631  2.438  2^M  2.26ft  >.186 
3.491  3.301  3.132  2.981  2.644  2.721  2.608  2.606  2.410  S.a»  2.244 
3.643  3.486  8.282  3.089  2.946  2.811  2.891  2.582  2.48ft  2.391  2  .*57 
3.811  3.666  3.386  3.211  3.061  2.910  2.782  £.666  2.660  2.446  2.874 


4.294  4.006  3.756  3.636 
4.669  4.2*  3.049  3.706 
4.876  4.498  4.177  3.904 
0.267  4.818  4.446  4.131 
6.797  5.216  4.771  4.403 


6.729  6.177  4.734 
6,427  6.706  6.151 
7.469  6.435  6.696 
9.223  7.632  6.460 
13.410  9.486  7.638 


14.668  9.830 

16.729 


3.344  3.173 
3.496  3.307 
3.666  3.469 
3.864  3.690 
4.094  3.880 


3.019  2.881 
3.129  2.990 
3^75  3.112 
8.427  3.249 
3.601  3.403 


3.306  3.679 
4.042  3278ft 
4.329  4.027 
4.684  4.321 
5.142  4.666 


2.757  2.644  ZJHX 
2.866  2.736  2.626 
2.966  2.636  2.717 
3.089  2.946  2.818 
S.S8  3.069  2.980 


2.060  1.969  1.932 

2.10e  2.039  1.979 1  1.90.  i<2A 

2.158  2.091  2.029 

2.217  2.147  2.082 


2.360  2.980 
£.433  2 .349 
2.511  2.421 
2.897  2.600 
2.690  2.687 


&jm.  3.209  3.065 

3.560  3.366  3.196 

3.766  3.647  3.366 

4J313  3.761  3.640 

4.SI9  4.015  3.797 


2.793  2.681 
2.907  2.786 
3.035  2.901 
3.180  3.030 
3.348  3.178 


5.181 
5.815 
8.794 
SB  8.519  |  6.949 
12.994  8.915 


*.422  7.377  6.222  5.407  4.901 

18.474  10.104  7.669  6.388  5.564 

28J8BB  11.044  9.089  6.S0Q 

12.431  8.612 

14.720 


2.406  |  2.387  E.8C 


2.679  2.487 

2.675  2.975 

2.700  2.872 

2.898  2.779 

3.C8C  2.898  |  2.780  3.53 


3.180  3.033 
3.365  3.167 
3.667  3.366 
3.000  3.574 
4.096  3.826 


4.475  4.136 
4.976  4.534 
5.694  5.069 
6.836  6.846 
9.135  7.137 


19.614  9.973 

43.870 


TABLE  2 


M.*  (^7)^ 


cosec  a 


29  30° 


2.116  2.C60 

2.171  2.102 

2.290  2.156 
2.217 


1.989  1.832 

2.038  I  1.878 
2.081 
2.147  I  2.0& 


2.207  2.136 
2.271  2.186 
2.939  2,2 63 
2.415  2.331 


saigai 


3.033  2.902 
3.187  3.0B9 
3.366  3.197 
3.S74  3.360 
3.820  3.896 


4.534 
.5.068 
6.646 
9.135  I  7.137 


4; 136  3.867 
4.534  4.182 

5.068  4.602 

6.846  5.173 

6.031 


2.402  2.324 
2.463  2.400 
2.873  2.462 
2.971  2.873 
2.780  2.973 


2.918 
3.065  I  2. 
3.230  3. 


1.962  1.911  1.864 
2.013  1.989  1.910 
2.097  X£!£  1.969 
2.12S  2.066  t.011 
2.187  2.124  2.096 


2.263  2.187  2.126 
2.324  2.263  2.188 
£.401  2.3E6  2.266 
2.464  2.404  2.330 
2.676  2.486  2.406 


1.819  1.777  1.736  1.701  1.667 
1.96S  1.330  1.779  1.740  1.706 
1.910  1.8SS  -1,80!  1.768  1.746 
1.960  1.913  1.808  1.6Z7  1.786 
2.012  1.965  1.919  1.972  UK 


2.016  1.966  1.921  1.979 
2.072  2.021  1.973  ui 
2.133  2.078  2.027  1.90 
2.187  2.140  2.066  2.066 
2.266  2JB06  2.148  2.016 


2.276  2.216  2.168 
2.«4  2.286  2.286 
2.489  2.367  2.30e 
2.681  2.464  2.882 
2.664  2.649  2.471 


1.446  1.420  1.397 
1.4J6  U90  1.486 
1.607  1.461  1.487 
1.640  1.614  1.4e9 


1.634  |  1.604  1.876  1.647  1.622 

1.643  1.610  1.696  1.666 

1.710  1.978  1.647  1.616  1.591 

i.76i  i.7i7  i.eee  i.665  i.tzi 

1.784  1.766  1.726  1.684  1.666 


2.068  2.014 

2.121  2.073 

2.137 

tjm  2.206 

2.397  2.378 


6.264  5.464  4.e90  4.466  4.131  3.140  3.016 


13.162  8.934 

17.083 


6.632  5.021  4.566 
6.864  6.8(7  6.173 
9.496  7.272  6.104 
at. 180  1C.T71  7.608 


3.667  5.384 
3.863  3.641 
4.218  3.948 


3.424  3.263  3.108  2.971  2.BS3  2.746  2.664  2.669  2.482  2.422  2.368  2.800  2.248 


2.388  2.328 
2.473  2.411 
2.673  2.604 
2.684  2.807 


2.810  2.728 
2.962  2.854 
8.116  3.003 
8.808  3.176 
5.659  3.380 


66°  67° 


72*  73*  74° 


78"  I  77*  I  78* 


1.152 
uiei 
utu 
1.250  I  1.241 


1.280  1.272 
1.312  1.304 
1.343  1.336 
1.376  1.368 
1.406  i.m 


1.851  1.846 

1.902  i.eoe 

1.956  1.949 

2.013  2.006 

2.073  2.066 


1.841  1.840  '1.842  1.846 
1.892  1.891  1.893  1.S98 


1.892  1.891  1.893  1.S98 
1.946  1.945  1.948  1.963 
2.003  2.002  2.005  2.011 
2.062  2.062  2.064  2.071 


ESHSIESSiS 

EE3ES3 

EE2JHEE3 


•5UI 


1.963 

1.069. 

2.014 

2.064 

2.079 

2.122 

2.147 

2.198' 

2.219 

2.269 

2.296 

2.361 

2.378 

2.437 

2.466 

2.631 

2.660 

2.631 

2.662 

2.740 

1.077'|  1.076 


1.124 

lilOB  I 
1.214 


1.287  1.307 

1.841  1.364 

1.886  1.482 

1.433  1.480 


1.423 

1.426 

1.432 

1.440 

1.460 

1.465 

1.479 

1.462 

1.466 

1.473 

1.482 

1.494 

1.609 

1.62$ 

1.601 

1.607 

1.614 

1.625 

1.639 

1.866 

1.577 

1.641 

1.648 

1.667 

1.669 

1.664 

1.604 

1.627 

1.568 

1.691 

1.601 

1.615 

1.682 

1.668 

1.679 

2.612  2.604 
2.612  2.716 
2.720  2.833 
2.839  2.996 


TABLE  3 


Mi=  (g^Tr)2  cosec(a-w) 


41  42°  43° 


1.631  1.49S  1.483  1.432 
1.6»  1.500  1.486  1.434 
1.541  1.604  1.470  1.436 
1.548  1.510  1.474  1.441 
1.567  1.E17  1.480  1.445 


1.866  1.610  1.566  1.526 


1.774  1.720  1.670  1.622  1.578  1.536  1.496  1.469 

1.734  1.738  1.688  1.637  1.591  1.547  1.506  1.467 

1.817  1.769  1.704  1.863  1.606  1.660  1.517  1.477 

1.842  1.781  1.724  1.671  1.621  1.574  1.630  1.468 


U461 


1.806  1.746  l.e91  I.eaB  1.58S  1.543  1.500 

1.770  1.712  1.356  1.607  1.5S9  1.514 

1.736  I.eaO  1.626  1.576  1.529 

1.7C6  1.648  1.596  1.54e 

1.729  1.671  1.617  1.56S 


1.759  1.697  1.640  1.586 
1.791  1.726  1.685  1.603 
1.826  1.757  1.694  1.631 
1.866  1.79e  1.725  1.662 
1.905  1.831  1.759  1.693 


2.011 
2.071 
2.13)  !  2. 
2.217  1 2. 


1.797  1.726 
1.639  1.764 
1.635  1.805 
1.933  1.651 
1.996  1.902 


2-306  2.177  2.062  1.930 


1.225  1.209  1.194  1.179  1.165  1.152  1.139  1.127  1.115  1.104  1.053  1.083  1.074  1.064  1.055  1.046  1.037 

1.215  1.190  1.183  1.167  1.153  1.13e  1.12S  1.112  1.099  1.087  1.076  1.064  1.064  1.043  1.033  1.C23  1.014 

1.206  1.189  1.173  1.156  1.141  l.UK  1.112  1.098  1.005  1.072  1.060  1.048  1.036  1.025  1.014  1.003  .9920 

1.199  1.181  1.164  1.147  1.130  1.115  1.100  1.066  1.072  1.068  1.045  1.032  1.020  1.008  .9963  .9837  .9722 


1.192  1.173  1.156  1.138 
1.186  1.167  1.148  1.130 
1.181  1.161  1.142  1.123 
1.177  1.156  1.136  1.117 
:  1.174  i.163  1.132  1.112 


1.150  1.128  1.108 

1.147  1.125  1.104 

1.148  1.123  1.102 
1.146  1.122  1.100 
1.145  1.121  1.099 


1.121  1.106  1.069 
1.113  1.098  1.080 
1.105  1.088  1.071 
1.099  1.081  1.064 
1.093  1.075  1.066 


1.068  J  1.069  1.061 
1.064  1.045 
1.060  1.041 
1.057  1.037 
1.056  1.034 


1.074  l.oeo  1.046  1.032  1.018  1.006 

1.064  1.049  1.034  1.(20  1.006  .9915 

1.065  1.039  1.024  1.009  .993e  .9790 

1.047  1.030  1.014  .9984  .9829  .9681 

1.C39  1.022  1.006  .9891  .9735  .9675 


1.032  1.016  .  9978  .  9609  .  9644  .  9483 

1.C27  1.009  .9911  .9736  ,9666  .9397 

1.C22  1.003  .  9649  .  9669  .  9496  .9321 

1.018  .9932  .  9794  .  9612  .  9433  .9252 

1.014  .0943  .9749  .9661  .9374  .9192 


1.171  1.146  1.122  1.098  1.075  1.063  1.032 


.9907  .9710 


.9327  .9138 


18  1.174  1.148  1.123  1.099  1.075  1.052  1.030  1.009  .9881  .9678  .9479  ,9284  .9091 

17°  1.177  1.160  1.124  1.100  1.075  1.052  1.030  1.007  .9860  .9663  .9448  .9248  .9063 

16°  1.182  1.154  1.127  1.101  1.077  1.052<  1.029  1.007  .9846  .9632  .9424  .9219  .9018 

19°  1.186  1.158  1.130  1.104  1.078  1.054  1.030  1.006  .9840  .9622  .9407  .9196  .8991 


.9920  .9794  .9665 
.9778  .9644  .9511 
.9649  .9507  ,9367 
.9628  .  9382  .9236 
.9422  .9(268,  .9116 


.9322  .9162  .9006 
.9233.  .9070  .8906 
.9150  .8964  .  8613 
.9077  .  6904  .8732 
.9013  .6834  .8656 


.6953  .8771  .8586 
.8902  .8714  .e528 
.8867  .8666  .8475 
.8619  .6622  .8426 
.8787  .  8584  .8886 


.9642 
.9376 
.9225  ,9037 
.9068  .8910 
.8962  .  8810 


1.C21  1.014  1.007 

.9949  .9857  .9764 

.(709  .9603  .9493 
.9490  .9369  .9253 


.9160  .9031 
.8974  .8833 
.8600  .8864 
.8642  .8487 
.8490  .6338 


.8668  .8529 
.8677  .8412 
.8477  .0307 
.8384  .8211 
.8301  .8122 


.6225 
.8155 
.8094 
.8040  I  .7044 
.7990  .7791 


,8364  .  8200  1.803011*2 
.8245  .8074 
.8126  .7969 
.8034  .7053 
.7942  .7768 


.7066  .7669 
.7780  .7588 
.7711  .7514 
.7648  .7449 
.7691  .7388 


1.163  1.134  1.107 


1.068  1.031  1.007  .9940 


.8966  .6700  .6564  .8350  .8149  .7946  .7743  .7540  .7334  .712 


1.169  1.140  1.112  1,084  1.058  1.033  1.008  .9846  .9616  .9389 

1.176  1,145  1.116  1.086  1.061  1.036  1.010  .  9e61  .9622  .  9391  .9164 
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3.065  3.091  3.124  3.162  3.208 
3.252  3.252  3.279  3.311  3.353 
3.457  3.483  3.519  3.660  3.613 
3.679  3.706  3.741  3.787  3.840 


1.065  1.0e5  1.069  j  1.072 

1.133  1.136  1.141 

1.202  1.208  1.215 

1.275  1.283  1.292  I  1.302 


1.350  1.361  1.373  1.386 
1.429  1.442  1.456  1.172 
1.511  1.527  1.543  1.563 
1.S97  1.614  1,635  ,1.658 
1.686  1.708  1.730  1.755 


2.331  2.371  2.415  2.464 
2.461  2.505  2.524  2.608 
2.6C2  2.646  2.700  2.759 
2.750  2.798  2.e56  2.922 
2.908  2.963  3.023  3.096 


3.079  2.961  3.203  3.281 
3.262  3.323  3.396  3.480 
3.401  3.459  3.526  3.605 
3.674  g.74?  3.831  3.929 
3.907  3.983  4.075  4.183 


1.157  1.165  1.172 
1.240  1.249  1.263 
1.327  1.342  1.358 


1.929  1.873  2.C20 
2.048  2.096  2.151 
2.172  2.228  2.288 
2.302  2.363  2.432 
2.440  2. 508  2.585 


2.587  2.661  2.743 
2.741  2.625  2.919 
2.906  2.997  3.103 
3.064  3.182  3.299 
3.272  3.381  3.506 


3.472  3.594  3.736 
3.690  3.824  3.977 
3.801  3.926  4.240 
4.180  4.340  4.526 
4.457  4.632  4.839 


1.0S3  1.C9&  1. 10- 

1.191  1.233  3.216 

1.233  1.31!  1.331 

1.396  1.423  1.452 


1.503  1.540  1.5'1'7 

1.624  1.662  1.709 

1.743  1.7SB  1.648 

1.870  1.926  1.991 

2.0C2  2.068  2.143 


4.098  4.331  4.622 
4.380  4.642  4.965 
4.687  4.979  6.342 
5.020  5.349  5.753 
5.386  5.753  6.212 


9.422  10.465 


TABLE  5 


Functions  for  Calculations  of  Expansive.  Flow 


>* 

?w 

g(M) 

hW  j 

M 

f« 

K(m) 

1i(m) 

90° 

130.45° 

.5283 

.6339  j 

44° 

120.32° 

.2971. 

.4203 

89° 

130.45 

.5282 

.6339  ! 

43° 

119.54 

.2860 

.40S9 

88° 

130.45 

.5279 

.6336 

42° 

118.71 

.2746 

.3972 

87° 

130.45 

.5274 

.6332 

41° 

117.83 

.2630 

,3852 

86° 

130.45 

.5268 

.6326 

40° 

116.90 

.2512 

.3727 

85° 

130.44 

.5259 

.6319 

39° 

115.92 

.2391 

.3599 

84° 

130.44 

.5249 

.6310 

38° 

114. SS 

.2269 

.3467 

83° 

130..42 

.5236 

.6300 

37° 

113.77 

.2146 

.3331 

82° 

130.41 

.5222 

.6287 

36° 

112.61 

.2022 

.3192 

81° 

130.39 

.5206 

.6274 

35° 

111.37 

.1897 

.3050 

SOc 

130.37 

.5188 

.6258 

34° 

110.07 

.1772 

.2905 

79° 

*  130.34 

.5168 

.6241 

33° 

10S.70 

.1647  ‘ 

.2757 

78° 

130.31 

.5146 

.6222 

32° 

107.23 

.1522 

.2607 

77° 

-  130.27 

.5122 

.6201 

31° 

105.69 

.1399 

.2454 

76° 

130.22 

.5096 

.6178 

30° 

104. OS 

.127S 

.2301 

75° 

130.16 

.5067 

.6154 

29° 

102.36 

.1159 

.2145 

74° 

130.10 

.5037 

.6127 

28° 

100.56 

.1043 

.1990 

73° 

130.03 

.5005 

.6099 

27° 

9S.65 

.09313 

.1835 

72° 

129.95 

.4970 

.6069 

26° 

96.64 

.08237 

.1681 

71° 

129.85 

.4933 

.6037 

25° 

94.53 

.07211 

.1529 

70° 

129.75 

.4894 

.6003 

24° 

92.30 

.06243 

.1379 

69° 

129.64 

.JS53 

.5966 

23° 

89.96 

.0533S 

.12:33 

68° 

129.51 

.4809 

.5928 

22° 

S7.50 

.04502 

.1092 

67° 

i29.37 

.4763 

.5SS7 

21° 

84.91 

.03739 

.09502 

66° 

129.22 

.4715 

.5845 

20° 

82.19 

.03051 

.08274 

65° 

129.05. 

.4664 

.5S00 

19° 

79.34 

.02447 

.07063 

64° 

128.87 

.4611 

.5752 

18° 

76.34 

.01919 

.05937 

63° 

128.67 

.4555 

.5702 

.17° 

73.21 

.01169 

.04903 

62° 

128.45 

.4497 

.5650 

16° 

69.94 

.01095 

.03977 

61° 

128.22 

.4436 

.5595 

15° 

66.51 

.00791 

.03153 

60° 

127.97 

.4372 

.5538 

14° 

62.95 

.00552 

.02438 

59° 

127.69 

.4306 

.5478 

13° 

59.23 

.00370 

.01832 

58° 

127.40 

.4237 

.5415 

12° 

55.37 

.00237 

-01331 

57° 

127.08 

.4165 

.5349 

11° 

51.37 

.00143 

.00931 

56° 

126.74 

.4090 

.5281 

10° 

47.22 

.00081 

.00621 

55° 

126.38 

.4013 

.5209 

9° 

42.94 

.00043 

.00392 

54° 

125.99 

.3933 

.5135 

8° 

38.53 

.00020 

.00232 

53° 

125.57 

.3850 

.5057 

7° 

34.00 

.00009 

.00126 

125.12 

.3763 

.4976 

6° 

29.36 

.00003 

.00061 

51° 

124.64 

.3675 

.4891 

5° 

24.63 

.00001 

.00026 

50° 

124.14 

.3583 

.4804 

4° 

19.81 

.00000 

.00009 

49° 

123.59 

.3488 

.4713 

3° 

14.92 

.00000 

.00002 

48° 

123.02 

.3390 

.4618 

2° 

9.98 

.00000 

.00000 

47° 

122.40 

.3289 

.4520 

i° 

5.00 

.00000 

.00000 

46° 

45° 

121.75  ' 

121.06 

15 

.3186 

.3080 

.4418 

.4312 

0° 

0.00 

.00000 

.00000 

TABLE  6 

Wedge  Angle  for  Sonic  Flow  Behind  Shock  Wave  («,).  Maximum  Wedge  Angle  (uw) 


M. 

w. 

w» 

Mo 

w. 

Win 

1.01 

0.05° 

1.75 

17.70° 

IS. 12° 

1.02 

0.13 

1.76 

17.DS 

IS.  34 

1.03 

0.24 

•  1.77 

IS. 20 

15:55 

1.04 

0.37 

If-  t  Jr.?5® 

1.78 

18.41 

18.76 

1.05 

0.52 

«fia»  » 

1.79 

18.63 

18.97 

1.G6 

0.67 

0.73 

1.80 

IS. 84 

19.18 

1.07 

0.S4 

0.91 

1.81 

19.05 

19.39 

1.0S 

1.02 

1.10 

1.82 

19.26 

19.59 

1.09 

1.21 

1.30  ! 

1.83 

19.46 

19. SO 

1.10 

1.41 

1.51  ' 

1.S4 

19.67 

20.00 

1.11 

1.61 

1.73  1 

1.S5 

19.87 

20.20 

1.12 

1.82 

1.96 

1.86 

20.07 

20.40 

1.13 

9.  04 

2.19 

1.87 

20.27 

20.59 

1.14 

2.26 

2.43 

l.SS 

20.47 

20.78 

1.15 

2.49 

2.67 

1.89 

20.67 

20.98 

1.16 

2.73 

2.92 

1.90 

.20.80 

21.17 

1.17 

2.97  1 

3.17 

1.91 

21 .05 

21,30 

1.18 

3.21 

3.42 

1.92 

21.24 

21.54 

1.19 

3.45 

3. 6S 

1 .93 

21 .43 

21.73 

1.20 

3.70 

3.94 

1.61 

21.62 

21.91 

1.21 

3.95 

4.21 

1 .95 

21 ,8l 

22.09 

1.22 

4.21 

4.48 

1.96 

21.99 

22.27 

1.23 

4.40 

4.74 

1.97 

22.17 

22.45 

1.24 

4.72 

5.01 

1.98 

22.35 

22.63 

1.25 

4.99 

5.29 

1.99 

22.53 

22. SO 

1 .20 

5.25 

5.56’ 

2.00 

22.71 

22.97 

1.27 

5.52 

5. S3 

2.01 

22. S-S 

23.14 

1.28' 

5.78 

6.11 

2.02 

23.05 

23.31 

1.29 

6.05 

6.39 

2.03 

23.23 

23.48 

1.30 

6.32 

6.66 

2.04 

23.40 

23.65 

1.31 

6.59 

6.94 

2.05 

23.56 

23.81 

1.32 

6.86 

7.22 

,2.06 

23.73 

23.9S 

1.33 

7.13 

7.49 

2.07 

23.90 

24.14 

1.34 

7.40 

7.77 

2. OS 

24.06 

24.30 

1.35 

7.67 

8.05 

2.09 

24.22 

24.46 

1.36 

7.94 

8.32 

2.10 

24.38 

24.61 

1.37  • 

8.21 

8.60 

2.11 

24.54 

24.77 

1.38 

8.43 

8.88 

2.12 

24.  ■’O 

24.92 

1.39 

8.76 

9.15  ' 

2.13 

24.85 

25.08 

1.40 

9.03 

9.43 

2.14 

25.01 

25.23 

1.41 

9.30 

9.70 

2.15 

25.16 

25.3S 

1.42 

9.57 

9.97 

2.16 

25.31 

25.52 

1.43 

9.84 

10.25 

2.17 

25.46 

25.67 

1.44 

10.10 

10.52 

2.18 

25.61 

25.82 

1.45 

10.37 

10.79 

2.19 

25.76 

25.% 

1.46 

10.64 

11.05 

2.20 

25>90 

26.10 

1.47 

10.90 

11.32 

2.21 

26.05 

26.24 

1.48 

11.17  ’ 

11.59 

2.22 

26.19 

26.38 

1.49 

11.43 

11.85 

2.23 

26.33 

26.52 

1.50 

11.69 

12.11 

2.24 

26.47 

26.66 

1.51 

11.96 

12.37 

2.25 

26.61 

26.80 

1.52 

12.21 

12.63 

1.53 

12.47 

12.89 

2.3 

27.28 

27.45 

1.54 

12.73 

13.15 

2.4 

28.53 

28.68 

1.55 

12.99 

13.40 

2.5 

29.67 

29.80 

1.56 

13.24 

13.66 

2.6 

30.70 

30.81 

1.57 

13.49 

13.91 

2.7 

31.64 

31.74 

1.58 

13.74 

14.16 

2.8 

32.50 

32.59 

1.59 

13.99 

14.41 

2.9 

33.29 

33.36 

1.60 

14.24 

14.05 

3.0 

34.01 

34.07 

1.61 

14.49 

14.90 

1.62 

14.73 

15.14 

,  4 

3S.75 

38.77 

1.63 

14.98 

15.38 

5 

41.11 

41.12 

1.64 

15.22 

15.62 

6 

42.44 

42.44 

1.65 

15.46 

15.86 

1.66 

15.70 

16.09 

8 

43.79 

43.79 

1.67 

15.93 

16.32 

1.68 

16.17 

16.55 

10 

44.43 

44.43 

1.69 

16.10 

1G.7S 

. 

1.70 

16.63 

17.01 

15 

45.07 

45.07 

1.71 

16.86 

17.24 

1.72 

17.09 

17.46 

20 

45.29 

45.29 

1.73 

17.31 

17.68 

1.74 

17.54 

17.90 

CO 

♦ 

45.58 

45.58 

im if. 


